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Abstract. The description of baryons as fully relativistic bound states of quark and glue reduces to an 
effective Bethe-Salpeter equation with quark-exchange interaction when irreducible 3-quark interactions are 

os : 

neglected and separable 2-quark (diquark) correlations are assumed. This covariant quark-diquark model 
of baryons is studied with the inclusion of the quark substructure of the diquark correlations. In order to 
maintain electromagnetic current conservation it is then necessary to go beyond the impulse approximation. 
A conserved current is obtained by including the coupling of the photon to the exchanged quark and direct 
"seagull" couplings to the diquark structure. Adopting a simple dynamical model of constituent quarks and 
exploring various parametrisations of scalar diquark correlations, the nucleon Bethe-Salpeter equation is 
solved and the proton and neutron electromagnetic form factors are calculated numerically. The resulting 
magnetic moments are still about 50% too small, the improvements necessary to remedy this are discussed. 



The results obtained in this framework provide an excellent description of the electric form factors (and 
charge radii) of the proton, up to a photon momentum transfer of 3.5GeV 2 , and the neutron. 



PACS. 11. 10. St (Bound states; Bethe-Salpeter equations) - 13.40.Gp (Electromagnetic form factors) 
14.20.Dh (Protons and neutrons) 
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1 Introduction 

To the precision accessible to current measurements, the 
proton is the only known hadron that is stable under the 
effect of all interactions. Protons are thus ideal for use in 
beams or targets for scattering experiments designed to 
explore the fundamental dynamics of the strong interac- 
tion. From this, it follows that more is known about the 
proton and its nearly stable isospin partner, the neutron, 
than any other hadron. There is an abundance of observ- 
able properties of the nucleon, from elastic scattering form 
factors to electromagnetic form factors and parton distri- 
butions, which are being measured with increasing pre- 
cision at various accelerator facilities around the world. 
In particular, the ratio of the electric and the magnetic 
form factor of the proton is subject to current measure- 
ments at TJNAF, and current experiments at MAMI |jj 
and NIKHEF show great promise that soon we will 
have a precise knowledge of the electromagnetic proper- 
ties of neutron as well. Hence, the development of an ac- 
curate and tractable covariant framework for the nucleon 
in terms of the underlying quarks and gluons is clearly 
desirable. 

Models of the nucleon and other baryons are numerous 
and have had varying degrees of success as they are usu- 
ally designed to describe particular properties of baryons. 
Some of the frameworks that have been employed are 
non-relativistic f|,(|,D an< ^ re l a tivistic |J quark poten- 
tial models, bag models jj],D, skyrmion models P, 10|,pT|, 
|l2|| or the chiral soliton of the Nambu- Jona-Lasinio (NJL) 
model [^3[|lJ] . The relativistic bound state problem of 3- 
quark Fadd eev ty pe was stud ied extensively within the 

and its non-local general- 
In addition, some 



NJL model |15 



19,2' 



ization, the global color model |2l|,[22f 
complementary aspects of these models have been com- 
bined, e.g., the chiral bag model (2^,Q and a hybrid 
model that implements the NJL-soliton picture of baryons 
within the quark-diquark Bethe-Salpeter (BS) framework 
©■ 

The present study is concerned with the further de- 
velopment of a description of baryons (and the nucleon, 
in particular) as bound states of quarks and gluons in a 
fully-covariant quantum field theoretic framework based 
on the Dyson-Schwinger equations of QCD. Such a frame- 
work has already reached a high level of sophistication for 
mesons. In these studies, the quark-antiquark scattering 
kernel is modeled as a confined, non-perturbative gluon 
exchange. Such a gluon exchange can be provided by so- 
lutions to the Dyson-Schwinger equations for the propaga- 
tors of QCD in the covariant gauge p6| , |2^ | . In phenomeno- 
logical studies the gross features of such solutions are mim- 
icked by the use of a phenomenological quark-antiquark 
kernel. With this kernel, the dressed-quark propagator is 
obtained as the solution of its (quenched) Dyson-Schwinger 
equation (DSE), and the same kernel is then employed in 
the quark-antiquark Bethe-Salpeter equation (BSE)[^]from 
which one obtains the masses of the meson bound states 



1 We warn the reader not to confuse the meaning of this 
acronym with another one adopted recently. See, e.g., Ref. |2q|. 



and their BS amplitudes. Once the dressed-quark propa- 
gators and meson BS amplitudes have been obtained, ob- 
servables can be calculated in a straight-forward manner. 
The success of this approach has been demonstrated in 
many phenomenological applications, such as electromag- 
netic form factors p9[ , decay widths, it-it scattering j30| , 
vector-meson electroproduction [ pl[ , to name a few. Sum- 
maries of this Dyson-Schwinger/Bethe-Salpeter descrip- 
tion of mesons may be found in Refs. |3^,|33|,|3~4||. 

The significantly more complicated framework required 
for an analogous description of baryons based on the quan- 
tum field theoretic description of bound states in 3-quark 
correlations has meant that baryons have received far less 
attention than mesons. Consequently, the utility of such 
a description for the baryons is considerably less under- 
stood, even on the phenomenological level. (For example, 
the ramifications of various truncation schemes, required 
in any quantum field theoretic treatment of bound states, 
which have been explored extensively in the meson sector 
are completely unknown in the baryon sector.) The devel- 
opments in this direction employ a picture based on sepa- 
rable 2-quark correlations, i.e., diquarks, interacting with 
the 3rd quark which allows a treatment of the relativistic 
bound state in a manner analogous to the 2-body BS prob- 
lem |35[[56],|37]]. Although these studies employ simplified 
model assumptions for the quark propagators and diquark 
correlators, they demonstrate the utility of the approach 
in general, and these simplified model assumptions can be 
replaced by more realistic ones in the future as more is 
understood about the underlying dynamics of quarks and 
gluons. 

In the present article we generalize this framework for 
the 3-quark bound state problem in quantum field the- 
ory to accommodate the quark-substructure of the sepa- 
rable diquark correlations. This necessary extension has 
important implications on the electromagnetic properties 
of the nucleons. We explicitly construct a conserved cur- 
rent for their electromagnetic couplings, and we verify an- 
alytically that it yields the correct charges of both nucle- 
ons. To achieve this we employ Ward and Ward-Takahashi 
identities for the quark correlations which arise from elec- 
tromagnetic gauge invariance. In order to avoid unneces- 
sary complications, our present study considers only sim- 
ple constituent-quark and diquark propagators. However, 
the generalization to include dressed propagators, most 
importantly to account for confined nature of quarks and 
diquarks, requires only minor modifications. 

The organization of the article is as follows. In Sec. @, 
some general properties of the two-quark correlations (or 
diquarks) are discussed. The importance of constructing 
a diquark correlator that is antisymmetric under quark ex- 
change is emphasized. The parametrisations of the Lorentz- 
scalar isoscalar diquark correlations which are employed 
in the numerical calculations of the subsequent sections 
are introduced to reflect these properties in deliberately 
simple ways. The framework is sufficiently general to ac- 
commodate the results for the diquark correlators as they 
become available from studies of the underlying dynam- 
ics of quark and gluon correlations in QCD. Of course, 
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diquark correlations other than the scalar diquark will be 
important for a more complete description of the nucleoli. 
Certainly necessary is the inclusion of other channels, such 
as the axialvector diquark, when the description is ex- 
tended to the octet and the decuplet baryons. Nonethe- 
less, for the purposes of the present study, only scalar di- 
quarks are considered; the generalization to include other 
diquarks is straight-forward jj7|. In Sec. |3[ the nucleon BS 
equation and the quark-exchange kernel are introduced. 
One particular feature of this kernel is that it necessar- 
ily depends on the total momentum of the nucleon bound 
state P n . This conclusion is based on general arguments, 
such as the exchange-antisymmetry of the diquark corre- 
lations. It is important to obtain the correct normaliza- 
tions and charges of the nucleon bound states. In partic- 
ular, to ensure current conservation requires a consider- 
able extension beyond the impulse approximation in the 
calculations of electromagnetic form factors for baryons.^ 
The numerical solutions for the nucleon BS amplitudes 
obtained herein preserve the invariance of observables un- 
der a re-routing of the relative momentum, a requirement 
that follows from the translational invariance of the rela- 
tivistic bound state problem. In Sec. |], the normalization 
condition for the nucleon BS amplitude is derived and the 
nucl eon electromagnetic current is obtained in Sec. H In 



Sec. 5.1, the "seagull" contributions to the electromag- 
netic current, necessary for current conservation, are de- 
rived from the Ward and Ward-Takahashi identities. In 
Sec. ^, expressions for the electromagnetic form factors 
of the nucleon ar e de rived, the numerical calculations are 
described in Sec. and the results for the form factors 



are presented and discussed in Sec. 6.2. The conclusions 
of this study are provided in Sec. [7. Several appendices 
have also been included with further details which may 
provide the interested reader with addition insight into 
the framework. 



2 Diquark Correlations 

To obtain a solution of the 3-particle Faddeev equations in 
quantum field theory requires a truncation of the interac- 
tion kernel. A widely-employed truncation scheme is to ne- 
glect contributions to the Faddeev kernel which arise from 
irreducible 3-quark interactions. This allows one to rewrite 
the Dyson series for the 3-quark Green function as a cou- 
pled system of equations. The first being the BS equation 
for the 2-quark scattering amplitude and the other being 
the Faddeev equation, which describes the coupling of the 
3rd quark to these 2-quark correlations. As a result, the 
nature of the gluonic interactions of quarks enters only 
in the BSE of the 2-quark subsystem, i.e., in the quark- 
quark interaction kernel. The solution of the full inhomo- 
geneous BS problem for this 2-quark system is simplified 
by assuming separable contributions (explained below), 



The P„-dependence of the exchange-kernel violates one of 
the necessary conditions for current conservation in the impulse 
approximation. See Sec. |^. 



hereafter referred to as diquarks, to account for the rele- 
vant 2-quark correlations within the hadronic bound state. 
Herein, a "diquark correlation" thus refers to the use of a 
separable 4-quark function in the 3 representation of the 
SU(3) color group. The utility of such diquark correla- 
tions for a description of baryon bound states is a central 
element of the present approach. 

While in the NJL model the 2-quark scattering ampli- 
tude has the property of being separable, in general this 
is an additional assumption useful to simplify the 3-quark 
bound state problem in quantum field theory. An example 
for separable contributions would be (a finite sum of) iso- 
lated poles at timelike total momenta P 2 of the diquark. 
Such poles allow the use of homogeneous BSEs to obtain 
the respective amplitudes from the gluonic interaction ker- 
nel of the quarks, and thus to calculate these separable 
contributions to the 2-quark scattering amplitude. 

However, the validity of using a diquark correlator 
to parametrise the 2-quark correlations phenomenologi- 
cally, does not rely on the existence of asymptotic diquark 
states. Rather, the diquark correlator may be devoid of 
singularities for timelike momenta, which may be inter- 
preted as one possible realization of diquark confinement. 
In principle, one may appeal to models employing a gen- 
eral, separable diquark correlator which need not have any 
simple analytic structure, in which case no particle inter- 
pretation for the diquark would be possible. The imple- 
mentation of this model of confined diquark is straight- 
forward and does not introduce significant changes to the 
framework. The use of diquark correlations in this capac- 
ity is quite general and does not necessarily imply the 
existence diquarks, which have not been observed experi- 
mentally. 

The absence of asymptotic-diquark states may be ex- 
plained in a number of ways. Although, it has been ob- 
served that solutions of the BSE in ladder approximation 
yield asymptotic color-3 diquark states p8|, when terms 
beyond ladder approximation are maintained, the diquarks 
cease to be bound p9| ]. That is, the addition of terms 
beyond the ladder approximation to the BS kernel, in a 
way which preserves Goldstone's theorem at every order, 
has a minimal impact on solutions for the color-singlet 
meson channels. In contrast, such terms have a signifi- 
cant impact on the color-3 diquark channels. In Ref. fS9[ , 
it was demonstrated that these contributions to BS ker- 
nel beyond ladder approximation are predominantly re- 
pulsive in the color-3 diquark channel. It was further- 
more demonstrated in a simple confining quark model that 
the strength of these repulsive contributions suffices to 
move the diquark poles from the timelike P 2 -axis and far 
into the complex P 2 -plane. While the particular, confin- 
ing quark model is in conflict with locality, the same effect 
was later verified within the NJL model |l0|]. This sug- 
gests that this mechanism for diquark confinement, which 
eliminates the possibility of producing asymptotic diquark 
states, might hold independent of the particular realiza- 
tion of quark confinement. 

In a local quantum field theory, on the other hand, col- 
ored asymptotic states do exist, for the elementary fields 
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momentum space as G™p s (p, q, P) where 



x 3 - 



-P8 



X2 
Xa 



Pj -Ps 
Fig. 2.1. The diquark pole in the 4-quark Green function. 



as well as possible colored composites such as the diquarks, 
but not in the physical subspace of the indefinite metric 
space of covariant gauge theories. The analytic structure 
of correlation functions, the holomorphic envelope of ex- 
tended permuted tubes in coordinate space, is much the 
same in this description as in quantum field theory with a 
positive definite inner product (Hilbert) space. In particu- 
lar, 2-point correlations in momentum space are generally 
analytic functions in the cut complex P 2 -plane with the 
cut along the timelike real axis. Confinement is interpreted 
as the observation that both absorptive as well as anoma- 
lous thresholds in hadronic amplitudes are due only to 
other hadronic states [[0]. However, the implementation 
of this algebraic notion of confinement seems much harder 
to realize in phenomenological applications. 

For the present purposes of developing a general frame- 
work for the description of baryon bound states, the ques- 
tion as to whether one should or should not model the 
diquark correlations in terms of functions with or without 
singularities for total timelike momentum P 2 is irrelevant. 
However, we reiterate that the present framework is able 
to accommodate both of these descriptions of confinement 
with only straight-forward modifications. 

The goal of the present study is to assess the utility 
of describing baryons as bound states of quark and di- 
quark correlations, in a framework in which the diquark 
correlations are assumed to be separable. (The term sepa- 
rable refers to the property that a 4-point Green function 
G(p, q; P) be independent of the scalar qp, where q and 
p are the relative momenta of the two incoming and out- 
going particles, respectively, and P is the total momen- 
tum.) To provide a simple demonstration of the general 
formalism, and its application to the calculation of nu- 
cleon form factors, we assume that the diquark correlator 
corresponds to a single scalar-diquark pole at P 2 = m 2 
which is both separable and constituent-like. In this de- 
scription, the (color-singlet) baryon thus is a bound state 
of a color-3 quark and a color-3 (scalar) diquark correla- 
tion. 

Assuming identical quarks, consider the 4-point quark 
Green function in coordinate space given by 



(2.1) 



(T(q 1 (x 3 )q a (x 1 )q (x 2 )qs{x4))) , 



where a, (3, 7, and S denote the Dirac indices of the quarks 
and T denotes time-ordering of the quark fields q a {x). The 
assumption of a separable diquark correlator, correspond- 
ing to the diagram shown in Fig. |^, can be written in 



z e 



APX 



G 



(xx,X2,X 3 ,X4) 



d{p) x~t a (p,P)xps(q,P) 



(2.2) 



where D(P) is the diquark propagator, X = ax\ + (1 — 
a)x3, Y = (l — a')x2 + a-'x^, the total momentum parti- 
tioning of the outgoing and incoming quark pairs are given 
respectively by er and a' both in [0, 1], and y = x\ — x%, 

Z = X2 — X4. 

As described above, the separable form of the diquark 
correlation of Eq. (^^) does not necessarily entail the ex- 
istence of an asymptotic diquark state. The framework 
developed herein makes no restrictions on the particu- 
lar choice of the diquark propagator D{P). Technically, 
model correlations which mimic confinement through the 
absence of timclike poles are easy to implement as shown 
in Ref. |§. 

Nonetheless, for the purpose of demonstration, here we 
employ the simplest form for this propagator correspond- 
ing to a simple (scalar) di quar k bound state in the 2-quark 
Green function G of Eq. (2.1). The appearance of such an 
asymptotic diquark state requires the diquark propagator 
to have a simple pole at P 2 — m 2 , where m s is the mass 
of the diquark bound state, i.e., 



D(P) = 



P 2 



(2.3) 



Then x(p>P) and its adjoint x(P>P) = 7oX t (P>- p bo are 
the BS wave functions of the (scalar) diquark bound state 
which are defined as the matrix elements of two quark 
fields or two antiquark fields between the bound state and 
the vacuum, respectively. Further details concerning those 
definitions are given in Appendix [A|. 

It is convenient to define the truncated diquark BS am- 
plitudes (sometimes referred to as BS vertex functions) x 
and x from the BS wave functions in Eq. (2.2) by ampu- 
tating the external quark propagators, 

X(p, P) ■= S-Hil-^P + p) X (P, P) S- lT (aP - p) , (2.4) 

X (p,P) := S- lT (aP-p) X (P,P) S- 1 ((l-a)P + p).(2.5) 

The convention employed here is to use the same symbols 
for both the BS wave functions and the truncated BS am- 
plitudes, with the latter denoted by a tilde. An important 
observation made in Appendix |X], which is of use in the 
following discussions, is that the for identical quarks, an- 
tisymmetry under quark exchange constrains the diquark 
BS amplitudes x to satisfy: 



x(p,p) = -x t (-p,p)\ 



(2.6) 



Note that a and (1 — a) have been interchanged here. 

In principle, at this point one could go ahead and spec- 
ify the form of the kernel for the quark-quark BSE in lad- 
der approximation, and obtain diquark BS amplitudes in 
much the same manner in which solutions for the mesons 
are obtained. However, as discussed above, the appear- 
ance of stable bound state solutions might be an artifact 
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of the ladder approximation rather than the true nature of 
the quark-quark scattering amplitude. Therefore, for our 
present purposes various simple model parametrisations 
for diquark BS amplitudes are explored, rather than us- 
ing a particular solution of the diquark homogeneous BSE. 
The motivation is to explore the general aspects and im- 
plications of using a separable diquark correlation for the 
description of the nucleon bound state. 



For separable contributions of the pole type (2.3), but 



possibly complex mass (with Re(m*) > 0), one readily 
obtains standard BS normalization conditions to fix the 
overall strength of the quark-quark coupling to the di- 
quark for a given parametrisation of the diquark struc- 
ture. These are obtained from the inhomogeneous quark- 
quark BSE for the Green function of Eq. (2.1) employing 
pole dominance for P 2 sufficiently close to rri 2 s R To sketch 
their derivation consider the inhomogeneous BSE which is 
of the general form, 



G{p,q,P) = (gM 1 (p,q,P)-K(p,q,P) 



(2.7) 



Here G^ denotes the antisymmetric Green function for 
the disconnected propagation of two identical quarks. The 
definition of its inverse G*(°) and a brief discussion of 
how to derive it, may be found in Appendix |X| With the 
simplifying assumption that the quark-quark interaction 
kernel K be independent on the total diquark momentum 
P; that is, K(p,q,P) = K(p,q) (which is satisfied in the 
ladder approximation for example), the derivative of G 
with respect to the total momentum P, gives the relation 



-^G(p, g> P)= / 



d*k d 4 k' 
(fcr) 5 (2^0* 



(2.8) 



G(p,k,P)[P» 



G<°> \k,k',P)\G(k',q,P). 



dP>* 

Upo n su bstitution of G(p, q, P) as given by Eqs. ( |2.2| ) 
and ( |2.3| ), and equating the residues of the most singu- 
lar terms, one obtains the non-linear constraint for the 
normalization of the diquark BS amplitudes \ an d X, 

d A p 



1 



Am 2 J (2tt) 4 
tr (s T ( P p)X(p,P) (p^-S( Pa ) ) x(p,P) 



(2.9) 



dP 



d 



+ tr ^x(p,P)S( Pa )x(p,P) [PfpS T ( Pp ) 

where p Q = p+(l— a)P andpp = —p+aP, i.e., P — p a +Pp 
and p = ap a — (l—a)pp. The scalar diquark contribution, 
relevant to the present study of the nucleon bound state, 
is color-3 and isosinglet. Lorentz covariance requires its 
Dirac structure to be the sum of four independent con- 
tributions, each proportional to a function of two inde- 
pendent momenta. For simplicity only a single term is 

3 For other separable contributions, e.g., of the form of non- 
trivial entire functions for which one necessarily has a singular- 
ity at \P 2 \ — ► oo , the use of the inhomogeneous BSE to derive 
normalization conditions relies on the existence of full solutions 



maintained here, which has the following structure: 

x(p,P) = J 5 C^P(p 2 ,pP) , (2.10) 

where C is the charge conjugation matrix (C = ij2^/o m 
the standard representation). The normalization constant 
N s is fixed from the condition given in Eq. ( [2.9] ) for a given 
form of P{p 2 ,pP). 

It may seem reasonable to neglect the dependence of 
this invariant function on the scalar pP; such a simplifi- 
cation would yield the leading moment of an expansion of 
the angular dependence in terms of orthogonal polynomi- 
als, which has been shown to provide the dominant con- 
tribution to the meson bound state amplitudes in many 
circumstances. However, in the present case the antisym- 
metry of the diquark wave function, c.f., Eqs. ( |2.6| ), entails 
that 



P(p 2 ,pP) = P(p 2 ,-pP) 



(2.11) 

7^ p, it is not 



>(1-<t) 

For a 7^ 1/2 and thus for p := p| ff ,j_ 

possible to neglect the pP dependence in the amplitude 
without violating the quark-exchange antisymmetry To 
maintain the correct quark-exchange antisymmetry, we as- 
sume instead that the amplitude depends on both scalars, 
p 2 and pP in a specific way. In particular, we assume the 
diquark BS amplitude is given by a function that depends 
on the scalar 



a) m s 

-2 



(1 - 2(j)pP -p 2 



(2.12) 



X := PaPfj - <x(l 

= (1 - 2a)pP-p- 

with p = (1 — a)p a — crpfj and P{ a ,0} as given above. For 
equal momentum partitioning between the quarks in the 
diquark correlation a = 1/2, the scalar x reduces to the 
negative square of the relative momentum, x = — p 2 . The 
two scalars that may be constructed from the available 



momenta p and P (noting that P 



is fixed) which 



have definite symmetries under quark exchange are given 
by the two independent combinations p a Pp (which is es- 
sentially the same as above x) and p 2 a — p\. The latter 
may only appear in odd powers which are associated with 
higher moments of the BS amplitude. Hence, these are 
neglected by setting 

P(p 2 ,pP) = P(x) (2.13) 

which can be s hown to satisfy the antisymmetry constraint 
given by Eq. fl2.1l| ) Vtr € [0,1]. Finally, the diquark BS 
normalization N s , as obtained from Eq. (2.9), is given by 



n: 



4 m? 



(27T) 4 

d 



P 2 {x) 



(2.14) 



P 7ns tr [S(P + (l-a)P)S(-p + crP)] 



dP 

The numerical results presented in the following sec- 
tions explore the ramifications of several Ansatze for P(x), 



-fn-p(^) 



7* 



or F EXP (a;) = exp {-ie/7exp}, (2.15) 



to (2.7) of the separable type. 



where the integer n — 1, 2... and corresponds to monopole, 
dipole,... diquark BS amplitudes. Their widths 7„, 7 EX p 
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are determined from the nucleon BSE by var ying them 
until the diquark normalization given by Eq. (2.14) and 
coupling strength g~ necessary to produce the correct nu- 
cleon bound-state mass are equal. This is carried out nu- 
merically and described in detail in Sec. ||. 

For completeness, various Gaussian forms that peak 
for values of x — xq > 0, 

exp{-(x-x ) 2 /7Lu 



}■ 



(2.16) 



are also explored. Such forms with finite xq were suggested 
for diquark amplitudes as a result of a variational calcu- 
lation of an approximate diquark BSE in Ref. and 
have been used in the nucleon calculations of Ref. ]22|. 
Therein, a fit t o the Gaussian form given by Ref. |}42| 
(and Eq. ( |2.16| ) above) was employed with a width of 
7 GAU ~ O.llGeV 2 and xo/7gau — F rom a calculation 
within the present framework, which is described in Sec. ||, 
we observe that the necessary value for 7 G au to obtain a 
reasonable nucleon mass is about an order of magnitude 
smaller than the value given in Ref. p^ j. Furthermore in 
Sec. [| we observe that the effect of a finite xq on the 
electric form factor of the neutron rules out the use of a 
Gaussian form with xo =/= 0. 

Some of the parametrisations used for diquark correla- 
tions in previous quark-diquark model studies of the nu- 
cleon |55[|3(],|37]], correspond to neglecting the substruc- 
ture of diquarks entirely; this may be reproduced in the 
present framework by setting P(p 2 ,pP) = 1. By neglect- 
ing the diquark substructure, the diquark BS normaliza- 
tions such as N s here are not well-defined. The strengths of 
the quark-diquark couplings are undetermined and chosen 
as free adjustable parameters of the model (one for each 
diquark channel maintained). In the present, more general 
approach, these couplings are determined by the normal- 
izations of the respective diquark amplitudes. At present 
this is the scalar diquark normalization N s alone which in 
turn determines the coupling strength g s = 1/N S between 
the quark and diquark which binds the nucleon. 

The aim of our present study is to generalize the notion 
of diquark correlations by going beyond the use of a point- 
like diquark and include a diquark substructure in a form 
similar to that of mesons. In the most general calculation 
of the three-body bound state, the strength of the quark- 
diquark coupling is not at one's disposal; it arises from 
the elementary interactions between the quarks. By us - 
ing the diquark BS normalization condition of Eq. (2.14), 
we can assess whether the quark-diquark picture of the 
nucleon is still able to provide a reasonable description 
of the nucleon bound state if the coupling strength is ob- 
tained from the diquark BSE rather than forcing the quark 
and diquark to bind by adjusting their couplings freely. 
Whether the quark-diquark coupling is sufficiently strong 
to produce bound baryons will eventually be determined 
by the strength of the quark-quark interaction kernel that 
leads to the diquark BS amplitude. 

Finally, the use of a diquark BS amplitude with a fi- 
nite width improves on the previous calculations in yet 
another respect. Without the substructure of the diquark, 
an additional ultraviolet regularization had to be intro- 
duced in the exchange kernel of the BSE for the nucleon 



in Refs. p5|,|36j ,|37j]. In the present study, the finite-sized 
substructure of the diquark leads to a nucleon BSE which 
is completely regular in the ultraviolet in a natural way. 

In this section, we have discussed some general fea- 
tures of the diquark amplitude employed in our present 
study. The important observations are the implications of 
its antisymmetry under quark exchange which constrains 
the functional dependence on the quark momenta, a nd th e 
derivation of the BS normalization condition, Eq. ( 2.14 ), 
to fix the quark-diquark coupling strength. By taking these 
into account, we will find that the precise form of the 
Ansatz for the diquark BS amplitude P(x) has little qual- 
itative influence on the resulting nucleon amplitudes as 
long as P{x) falls off by at least one power of x for large 
x corresponding to a large spacelike relative momentum. 



3 The Quark-Diquark Bethe-Salpeter 
Equation of the Nucleon 



By neglecting irreducible 3-quark interactions in the ker- 
nel of the Faddeev equation giving the 6-point Green func- 
tion that describes the fully-interacting propagation of 3 
quarks, the Dyson series for this 6-point Green function 
reduces to a coupled set of two-body Bethe-Salpeter equa- 
tions, see for example, Ref. [Q. As discussed in the previ- 
ous section, for the purpose of demonstrating the frame- 
work developed herein, we choose to employ constituent 
quark and diquark propagators. However, the framework 
itself is much more general. The assumptions under which 
it is developed require only that the irreducible 3-quark 
interactions are neglected in the nucleon Faddeev kernel 
and that the diquark correlations be well-parametrised by 
a sum of separable terms (which may or may not have 
poles associated with asymptotic diquark states). 

Maintaining only the (flavor-singlet, color-3) scalar di- 
quark cha nnel in these correlations, corresponding to Eqs. 
(2.2) and (2.3), the appearance of a stable nucleon bound 
state coincides with the development of a pole in the Green 
function G a p describing the fully-interacting (spin-1/2) 
quark and (scalar) diquark correlations which is of the 
form, 



G P S(P>k,P n ) 



(3.1) 



PS - Ml + 



-rl>(-k,P n ) 



a/3 



Here k and p are the relative momenta between the quark 
and diquark, incoming and outgoing, respectively, P n is 
the four- momentum of the nucleon with mass M n , a and (3 
denote the Dirac indices of the quark, and the nucleon BS 
wave function ip(p, P n ) is related to it's adjoint according 
to 

$(p, P n ) = Tbl^C-P. Pnho = C^ T {p, -P^C- 1 . (3.2) 
The truncated nucleon BS amplitude tp is defined as 
^(p,P„) - D((l- v )P n -p)S( V P n +p)$(p,P n ) , (3.3) 
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Fig. 3.1. The quark-exchange kernel of the nucleon BSE. 



where 77 S [0, 1] is the fraction of the nucleon momen- 
tum P n carried by the quark. The resulting homogeneous 
Bethe-Salpeter equation for the nucleon bound state reads, 

~ f d 4 k 

i> a p(p, Pn) = J ^— ^ K ai (p, k; P n ) ip^piji, P n ) , (3.4) 

where K a p(p, fc; P n ) is the kernel of the nucleon BSE which 
represents the exchange of one of the quarks within the 
diquark with the spectator quark. Maintaining the quark- 
exchange antisymmetry of the diquark correlations, this 
kernel provides the full antisymmetry within the nucleon 
(i.e., Pauli's principle) in the quark-diquark model ]15fl . 
The exchange kernel is shown in Fig. [| and given by 

K a0 (p,k;P n ) = (3.5) 

1 ~T 

- 2 XayiPl, q + Pa) S^ s (q) X S0 (p 2 ,q + Pfi) 

= ^ P(x 1 )P(x 2 ) S a0 (q) , (3.6) 

where the factor -1/2 arises from the flavor coupling be- 
tween the quark and diquark and p\ and p 2 are the relative 
momenta of the quarks within the incoming and outgoing 
diquark, respectively, such that 

Pi = cp a - (l-o - )? , and p 2 = (l-a')q - a'pp . (3.7) 

The respective momentum partitionings are a and a' and 
need not be equal. The total momenta of the incoming and 
the outgoing diquark are q+p a and q+pp- The momentum 
q of the exchanged quark is expressed in terms of the total 
nucleon momentum P n and relative momenta k and p as 



(1-27?)P» 



P- 



(3.8) 



Using the definitions of q and the quark momenta p a — 
r\P n + p and pp — r\P n + k, the relative momenta in the 
diquark BS amplitudes can be expressed as 

pi = {ar]-{l-(T)(l-2r]))P n + (l-a)k + p, (3.9) 
P2 = {{l-a'){l-2r,)-a'ri)P n -k-{l-a')p. (3.10) 

The corresponding arguments of the quark-exchange-anti- 
symmetric diquar k BS amplitudes follow readily from their 
definition in Eq. (2.12), 

xx = -p\ - {l-2a){{l-rj)p x P n - p x k) , (3.11) 
(l-2a')((l-r ] )p 2 P n -p 2 p) . (3.12) 



•7'2 



-Pi 



Before proceeding, it is worth summarizing some of the 
important aspects of this framework: 



1. The dependence of the exchange kernel on the total 
nucleon bound-state momentum P n is crucial in order 
to obtain the correct relationship between the elec- 
tromagnetic charges of the proton and neutron bound 
states and the normalizations of their BS amplitudes. 

2. The momentum q of the exchanged quark is indepen- 
dent of the nucleon momentum P n only for the partic- 
ular value of momentum partitioning rj = 1/2. 

3. The relative momenta p\ and p 2 of the quarks within 
the diquarks are only independent of the total momen- 
tum P n of the nucleon for the particular choice: 



1 - 2n 



1 



V 



(3.13) 



The symmetric arguments of x\ and x 2 in the diquark 
BS amplitudes are independent of the total momentum 
P n only if, in addition to the above criteria, a = a 1 = 
1/2 as well, and hence 77 = 1/3. In fact, this conclusion 
can be further generalized: Any exchange-symmetric 
argument of the diqua rk amplitude can differ from the 
definition of Eq. ( 2.12 ) only by a term that is propor- 
tional to the square of the total diquark momentum. 
From this, it is possible to show that the diquark BS 
amplitudes can be independent of P n only if 77 = 1/3. 
This is shown in Appendix [b|. 

Unlike the ladder-approximate kernels commonly em- 
ployed in phenomenological studies of the meson BSE, 
the quark-exchange kernels of the BSEs for baryon bound 
states, such as the one depicted in Fig. ||, necessarily de- 
pend on the total momentum of the baryon P n for all 
values of r\ £ [0,1]. This has important implications on 
the normalization of the bound-state BS amplitudes as 
well as on the calculation of the electromagnetic charges 
of the bound states. In particular, considerable extensions 
of the framework beyond the impulse approximation are 
required in the calculation of electromagnetic form fac- 
tors in order to ensure that the electromagnetic current is 
conserved for baryons. This issue is explored in detail in 
Sec. |. 

The form of the pole contribution arising from the nu- 
cleon bound stat e to the quark-diquark 4-point correlation 
function in Eq. (3.1) constrains the nucleon BS amplitude 
to obey the identities: 



4>( P ,Pn)A+(P n ) =^(p,P n ) , 

A+(P n )$(p,P n ) = ^(p,P„) , 



(3.14) 
(3.15) 



where A+(P n ) = (#„ + M n )/2M n . It follows from this 
that the most general Lorentz-covariant form of the nu- 
cleon BS amplitude can be parametrised by 

^(p,Pn) = (3.16) 

Si(p,P n )A+(P n ) + S 2 ( P ,P n )S(p,P n )A+(P n ) , 

$(p,P n ) = (3.17) 
Si(-p,P n )A+(P n ) + S 2 {- P ,P n )A+(P n )S(~p,P n ) , 

where S(p,P n ) = (jj> - pP n /M n )/M n , and Si(p,P n ) and 
S 2 (p, P n ) are Lorentz-invariant functions of p and P n . This 
provides a separation of the positive and negative energy 
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Fig. 3.2. The first 5 moments of the nucleon BS amplitudes 
5*1 (left) and S 2 (right). 



components of the nucleon BS amplitude. In Appendix |b], 
some of the consequences of this decomposition are ex- 
plored. In particul ar, it allows one to rewrite the homoge- 
neous BSE of Eq. (3.4) in a compact manner, in terms of 
a 2-vector S T (p,PnT= (Si(p,P„), S 2 (p,P„)) as 

S(p,P n )= (3.18) 
' 1 '' k P(x 1 )P(x 2 )D(k s )T{p,k,P n )S(k,P n ), 



2N$J (2tt)' 

where k s = (l—r/)P n — k and T(p, k, P n ) is a 2 x 2 matrix 
in which each of the four elements is given by a trace 
over the Dirac indices of the quark propagator S(k q ) (with 
k q = rjP n + k), the propagator S(q) of the exchanged 
quark and a particular combination of Dirac structures 
derived from the decomposition of the nucleon amplitude 
in Eq. ( 3.16 ), see Appendix [b|. 

Upon carrying out these traces, the resulting BSE is 
transformed into the Euclidean metric by introducing 4- 
dimensional polar coordinates corresponding to the rest 
frame of the nucleon according to the following prescrip- 
tions (where "— >" denotes the formal transition from the 
Minkowski to the Euclidean metric): 

P 2 - -P 2 , P 2 -> M 2 n , pP n - iM n py , (3.19) 
and S(p,P n ) -> S(p,y) . 

Written in terms of these variables, the nucleon BS ampli- 
tude is a function of the square of the relative momentum 
p 2 and the cosine y G [—1, +1] of the azimuthal angle be- 
tween the four- vectors p and P n . The dependence of the 
nucleon BS amplitudes S(p, y) on the angular variable y 
is approximated by an expansion to the order A in terms 
of the complete set of Chebyshev polynomials T n (y), 



N-l 



S(p,y) ~ S n (p)T n (y) 



71=0 



N 



Snip) = i n JJ^2 S (P' Vk)T n iVk) , 
fc=l 

V(fc- l/2) s 



(3.20) 
(3.21) 



where yt = cos 



A 



are the zeros of the Chebyshev polynomial T^iy) of de- 
gree N in y. Here, we employ Chebyshev polynomials of 
the first kind with a convenient (albeit non-standard) nor- 
malization for the zeroth Chebyshev moment from setting 





v 




J 







Fig. 3.3. Dependence of A/A^ s 2 on y for fixed a = a' = 1/2 
(solid) and with a, a' from Eq. (3.13) (dashed). The width 72 



is fixed to yield X/N s = 1 at 77 = 1/3. 

To = 1 / v/2. The explicit factor of (— i) n is introduced 
into Eq. ( [3.2C| ) to ensure that the moments S n {p) are real 
functions of positive p = \pp^ for all n. The BSE for these 
moments of the nucleon BS amplitude is then 

S m ip) = (3.22) 

k 3 dk 

i(p,k)S n (k) , 



1 



2N* J (An) 



N-l 

E 

71=0 



T 



where the (2 AT x 2 A) matrix T mn (p, k) is obtained from 
T(p, k, P n ) by expanding in terms of Chebyshev p olyno - 
mials both amplitudes S in the nucleon BSE of Eq. ( |3.18|) ; 
that is, summing over the yu on both si des o f Eq. (p.lSj) 



according to Eq. ( |3.2l| ) and using Eq. ( |3.20[ ) for the BS 
amplitude on the right-hand side. The definition of the 
matrix T mn appearing in Eq. (3.22) furthermore includes 



the explicit diquark propagator and BS am plitudes P(x\) 
and P{x 2 ) of the integrand in Eq. ( [3.18 ) as well as the 
integrations over all angular variables. Its exact form is 
provided in Appendix [Bj. 

In the calculations presented herein, we restrict our- 
selves to the use of free propagators for constituent quark 
and diquark, with masses m q and m s respectively, as the 
simplest model to parametrise the quark and diquark cor- 
relations within the nucleon. Measuring all momenta in 
units of the nucleon bound-state mass M n leaves m q /M n 
and m s /M n as the only free parameters. Using fo r the 
sca lar d iquark amplitude P(x) the forms of Eqs. ( [2.15| ) 
or ( 2.16| ) with fixed wid ths 7, the homogeneous BSE for 
the nucleon in Eq. ( |3.4|) is converted into an eigenvalue 
equation of the form, 



A^(>,Pn) 



d 4 k 



N 2 K(p,k,P n )ij{p,P n 



(3.23) 



with the additional constraint that A = N 2 (note that 
A^A" is independent of N s ). This equation is solved iter- 
atively and the largest eigenvalue A is found which cor- 
responds to the nucleon ground-state. To implement the 
cons traint , we calculate the diquark normalization A^ from 
Eq. (2.14) and compare it to the eigenvalue A. This pro- 
cedure is repeated with a new value for the width 7 of 
the diquark BS amplitude until the eigenvalue obtained 
from the BSE agrees with the normalization integral in 
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Fig. 3.4. The leading moments of the nucleon BS a mplitudes 
Si (left) and S2 (right) for the diquark amplitudes ( 2,15 2.16 ) 
with m, adjusted to yield Si,o(p)|j>=o.2M„ = 1/2. 



N 2 is satis- 



Eq. (|2.14D ; that is, until the constraint A 
fed. 

A ty pical solution of the nucleon BSE is shown in 
Fig. |3.2| . Five orders were retained in the Chebyshev ex- 
pansion. The figure depicts the relative importance of the 
first four Chebyshev moments of the nucleon amplitudes 
Si(p, Pn) an d Slip, P n ). The respective results for their 
fifth moments are too small (< 1 0~ 4 ) to be distinguished 
from zero on the scale of Fig. 3.2. This provides an indica- 



tion for the high accuracy obtained from the Chebyshev 
expansion to this order. This particular solution, which 
will be shown to provide a good description of the electric 
form factors in Sec. ^, was obtained for the d ipole form of 
the diquark amplitude, i.e., from Eq. (2.15) with n = 2, 
using m q — 0.62M n and m s = 0.7 M n (with a = a' = 1/2, 
T) = 1/3). The value for the corresponding diquark width 
72, necessary to yield A = N 2 (to an accuracy of 10~ 3 ), 
resulted thereby to be 72 = (0.294 M n ) 2 . 

The dependence of the BS eigenvalue on the momen- 
tum partitioning parameter 77 is shown in Fig. |3.3j . The 
complex domains of the constituent quark and diquark 
momentum variables, k 2 q and k 2 respectively, that are sam- 
pled by the integration over the relative momentum k in 
the nucleon BSE are devoid of poles for 1 — m s /M n < 77 < 
m q /M n . The pole in the mo men tum q 2 of the exchanged 
quark of the kernel in Eq. (|3.6| ) is avoided by imposing 
the further bound 77 > (1 — m q /M n )/2 on the momen- 
tum partitioning parameter. Therefore, with the present 
choice of constituent masses it is for values of 77 in the 
range 0.3 < rj < 0.6, for which our numerical methods can 
yield stable results. 

In principle, the integrations necessary to solve the nu- 
cleon BSE should always be real and never lead to an 
imaginary result (since m q + m s > M n ). More refined 
numerical methods would be necessary, however, when in- 
tegrations were to be performed in presence of the poles 
in the constituent propagators that occur within the inte- 
gration domain for values of 77 outside the above limits. 

The momentum partitioning within the nucleon 77 is 
not an observable. Hence, for any value of 77 for which 
our numerical results are stable, we expect the eigenvalue 
of the nucleon BSE A (which implicitly determ ines the 
nucleon mass) to be independent of 77. In Fig. |3.3| , the 
dependence on 77 of the nucleon BSE eigenvalue A is ex- 
plored using a fixed value for a = a' = 1/2 (solid curve) 



and using the values of a and a' from Eq. (3.12) (dashed 
curve). 

In the first case, the arguments of the diquark ampli- 
tudes simplify to Xi — —p 2 with pi = — (1 — 3r))P n /2 + 
p + k/2 and p 2 = (1 - drj)P n /2 - k - p/2. This implies 
that the real parts of Xi are guaranteed to be positive 
only for 77 = 1/3. For values of 77 ^ 1/3, a negative contri- 
bution arises from the timelike nucleon bound-state mo- 
mentum P n . If the n-pole forms for the diquark ampli- 
tudes are employed, this results in the appearance of ar- 
tificial poles whenever (1 — 3t7) 2 M^/4 > 7„. The value of 
72 = (0.294 M n ) 2 as it results here, would entail the ap- 
pearance of a pole for 77 > 0.53. This explains why our 
numerical procedure becomes unstable as 77 approaches 
the value 1/2 in this case. 

No such timelike contri but ion to the Xi arises in the 
second case shown in Fig. 3J3. Here, for a = a' = (1 — 
2?y) /( 1 — 77) the relative momenta pi are independent of 
P n , and t here are n o terms oc M 2 in xf, see, for exam- 
ple, Eqs. ( ^Tlpdl ). However, as 77 — > 1/2, we find that 
a, a' — > and the diquark normalization integrals are then 
affected by singularities. 

In the restricted range allowed to 77, the results for the 
nucleon BSE obtained herein are found to be independent 
of 77 to very good accuracy when the equal momentum 
partitioning between the quark in the diquarks is used, 
er = a' — 1 /2. In this case, the diquark normalization 
in Eq. ( 2.14 ) yields a fixed value for N 2 and any 77 de- 
pendence of the ratio X/N 2 has to arise entirely from the 
nucleon BSE. The 77 independence of A thus demonstrates 
that the solutions to the nucleon BSE are under good con- 
trol. 

The more considerable ^-dependence observed for a = 
a 1 = (1 — 277) / (1 — 77) arises from the dependence of N 2 on 
er, a'. Here, the limitations of the model assumptions for 
the diquark BS amplitudes are manifest. The dependence 
of the diquark BS amplitudes on the diquark bound-state 
mass, higher Chebyshev moments (which have been ne- 
glected) and the Lorentz structures (which we have ne- 
glected), are all responsible for this observed a depen- 
dence. Such a dependence would not occur, had the di- 



Fixed width Si(p)\ p = 
P{x) m a [M n ] 



0.2M„ = 

in a [M, 



1/2, see Fig 
] Vl [M n 



3.4 



] 9s = 1/N. 



n = l 
n = 2 
n — 4 
n — 6 
n = 8 
EXP 
GAU 



0.7 
0.7 
0.7 
0.7 
0.7 
0.7 
0.7 



0.685 
0.620 
0.605 
0.600 
0.598 
0.593 
(L572 



Fixed masses, see Fig 
n = l 0.7 
n = 2 0.7 
n = 4 0.7 
: 6 0.7 



n 
EXP 



0.7 



|3.5: 
TTB2 

0.62 
0.62 
0.62 
0.62 



0.162 
0.294 
0.458 
0.574 
0.671 
0.246 
0.238 

0.113 
0.294 
0.495 
0.637 
0.283 



117.1 
91.80 
85.47 
84.37 
83.76 
82.16 
71.47 

155.8 
91.80 
81.08 
78.61 
74.71 



Table 3.1. Summary of paramet ers used f or the various di- 
quark amplitudes P(x), c.f., Eqs. (2.15 2.16). 
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Fig. 3.5. The moments Si,o (left) and 52, o (right) of the nu 
cleon amplitudes for the diquark amplitudes of Eqs. ( 2.15 
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2.16|) 



quark amplitudes employed herein been calculated from 
the diquark BSE. 

For the numerical calculations of electromagnetic form 
factors presented in Sec. ||, we therefore choose to restrict 
the model to a = a' = 1/2 and vary ry. From the observed 
independence of the nucleon BSE solutions to r\ (when 
a = a' = 1/2), one expects to find that calculated nucleon 
observables, such as the electromagnetic form factors, will 
display a sim ilar independence of r\ as well. 

In Figure 3.4 the zeroth Chebyshev moments Sifiip) 
and 5*2,0 (p) obtained from a numerical solution of the nu- 
cleon BSE with a = a' = 1/2 and i] — 1/3 are show n, for 
various di quar k amplitudes of the forms given in Eqs. (|2~15|) 
and Eq. (2.16) with xq = 0. To provide a comparison be- 



tween these different forms of diquark BS amplitude, we 
have chosen to keep the diquark mass fixed and vary the 
quark mass until a solution of the nucleon BSE was found 
with the property that 



Si,o(p)|p=o.2Ar„ = I/ 2 



(3.24) 



The resulting values for m q along with the corresponding 
diquark wid ths and couplings g s = 1 /N s are summarized 



in Table 3.1 



We observe that the mass o f the quark required to sat- 
isfy the condition in Eq. ( |3.24| ) tends to smaller values for 
higher n-pole diquark amplitudes. On the other hand, for 
fixed constituent quark (and diquark) mass the higher n- 
pole diquark amplitudes lead to wider nucleon amplitudes. 
This is shown in Fig. 3.E and the corresponding diquark 
wid ths and couplings are given in the lower part of Ta- 
ble 3.1. The mass values here by c orrespond to the results 
shown in the previous Figure 3.4 for n = 2. 



The qualitative effec t of s hifting the maximum in the 
Gaussian forms in Eq. (2.16) for the d iqua rk amplitudes 
by a finite amount xo, is shown in Fig. 3.6. The curve for 
xq = resembles the Gaussian result given in Fig. |3.4| 



with masses m s = 0.7M n and m q = 0.572Af„, which are 
kept fixed in the results for finite shifts xq. The corre- 
sponding widths and normalizations of t hese forms for 
the diquark amplitudes are given in Table 3.2. While the 
widths 7 and the couplings g s are not free parameters 
in our approach, the additional parameter Xo/j GAU intro- 
duced into the form of the diquark BS amplitude is free 
to vary. In contrast to the Gaussian form of the model di- 
quark BS amplitude employed in Ref. |22]| with xo/"f GAU = 
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Fig. 3.6. Th e mo ment Si,o for the Gaussian diquark ampli- 
tudes of Eq. ( |2.16| ) with various shifts :to/7gau (left), see Ta- 
ble 3.2, and a comparison (right) of our result for ^o/tgau 
Mi II (with m s = 0.7M n , m q — 0.555A'/ n ) with a result we 
obtain from using 7gau = O.llGeV 2 as in Ref. |22|. 



0.19/0.11 ~ 1.73, by implementing the diquark normaliza- 
tion condition, we find that the width must be about 30 
times smaller to provide a sufficient interaction strength 
g s = necessary to bind the nucleon. 

To provide for a closer comparison with the results 
of Ref. p2|, we have also solved the nucleon BSE using 
the values of parameters employed in that study. That 
is, we have solved the nucleon BSE using P GA u(x) with 
the parameter xo/7gau = 0.19/0.11 and using the con- 
stituent quark and diquark masses m q = 0.555M n and 
m s — 0.7M n , respectively. (Note that we use the nucleon 
mass as the intrinsic momentum scale of our framework, 
whereas the momentum scale employed in Ref. |22] ] was 
1 GeV.) The value employed in Ref. ]2^] for the diquark 
mass m s = 0.568 GeV implies that M n = 811 MeV in 
order to compare to our calculations (with m s = 0.7 M n ), 
and our value for the quark mass (m q = 0.555M„) thus 
corresponds to m q = 0.45 GeV. With these parameters, 
we obtain 7gau = (0.216 M n f = 0.308 10" 
may be compared to 7 G au = (0.409 M n ) 2 
used in Ref. |E2l. 



2 GeV 2 , which 
= O.llGeV 2 



m s = 0.7M„, m q = 

GAU \/7GAU 



0.572M n : 

[M n ] g s 



1/N. 





1/4 

1/2 
3/4 

1 
5/4 
3/2 
7/4 



0.238 
0.209 
0.181 
0.154 
0.130 
0.109 
0.091 
0.077 



71.47 
69.09 
72.64 
81.96 
97.75 
121.3 
154.3 
198.5 



The same re- 
sults for the nu- 
cleon amplitu- 
des as shown 
above are used 
for the calcu- 
lations of the 
electromagnetic 
form factors of 
the proton and 
the neutron in 
Sec. |. We will 
show that when 
n-polc diquark 
BS amplitudes 
with n — 2 and n — 4 are employed, one obtains re- 
sults for the electromagnetic form factors that are in good 
agreement with the phcnomenological dipole fit for the 
electric proton form factor, while higher powers or expo- 
nential forms tend to produce form factors which decrease 
too fast with increasing momentum transfer Q 2 . We also 
show that use of a Gaussian form for the diquark BS am- 



Table 3.2. Width and normalizations of 
Gaussian amplitudes, 
ous shifts xq, see Fig. 



E. 



( Mq ) 

(right). 



for vari- 
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plitude with xq ^ leads to neutron electric form fac- 
tor that has a qualitatively different behavior than that 
obtained by using the other model forms of diquark BS 
amplitudes. Furthermore, this form leads to nodes in the 
neutron electric form factor for small Q 2 , a feature for 
which there is no experimental evidence. 



4 Normalization of Nucleon Bethe-Salpeter 
Amplitude 

In order to reproduce the correct electromagnetic charges 
of the physical asymptotic states, it is required that their 
Bethe-Salpeter (or Faddeev) amplitudes be normalized ac- 
cording to the normalization conditions obtained from the 
fully-interacting Green functions of the elementary fields. 
These normalization conditions are derived from inhomo- 
geneous BS (or Faddeev) equations. In the present frame- 
work, the quark-diquark 4-point Green function is given 
by (suppressing Dirac indices) 



G(p,q,P) = (GW \p,q,P)-K(p,q,P) 



(4.1) 



Here, K (p, q, P) is the quark- exchange kernel of Eq. ( p.6|) 
depicted in Fig. ||, and G^ (p, q, P) is the disconnected 
contribution to this Green function, given by 

G%{p,k,P) := (2^f5\p~k)S aP {k q )D{k s )5 m , , (4.2) 

with k q = rjP + k and k s = (l—rj)P — k. Taking the deriva- 
tive of G(p, k, P) with respect to the total momentum P 
and then examining the leading contributions in the limit 
that P 2 -> Ml which arise from the nucleon pole contri- 
bution from Eq. 



M n A+(P n ) = 



V 



one finds 

d 4 k 



(2tt) 4 (2tt) 4 

P,^rG-\ P ,k,p) 



(4.3) 



1>(k,P n ) 



p=p n 



The most important difference in normalizing the ampli- 
tudes in the Bethe-Salpeter framework developed here and 
a genuine two-body BSE, is the dependence of the BS ker- 
nel on the total momentum of the bound-state P. In phe- 
nomenological studies of 2-body bound states BS kernels, 
such as ladder-approximate ones for example, are most 
commonly employed in a form which does not depend on 
the bound-state momentum. However, in the present ap- 
proach, the original 3-body nature of the nucleon bound 
state requires the exchange kernel of the reduced BS equa- 
tion for quark and diquark to depend on the total momen- 
tum P n of the nucleon. In particular, when rj ^ 1/2, the 
propagator for the exchanged quark in the kernel depends 
on the total momentum P n and for rj ^ 1/3, the diquark 
BS amplitudes depend on P n . This added complication, 
which is easily avoided by ladder-approximate studies of 
meson bound states is unavoidable for studies of baryons. 
Thus, the normalization condition for the nucleon will al- 
ways contain contributions from the derivative of the ker- 
nel. 



Following the above procedure, one obtains a normal- 
ization condition of the form 

1 = (4.4) 

r}N q + (1 - V )N D + (1 - 2r])N x + (1 - S^Np , 
where 

d 4 k 



N q = 



pi* 



pn 



(2tt) 4 
tr 
d 4 k 



D(k s ) 



(4.5) 



d 



i>(-k,P) [ g^S(k q ))ip(k,p) 



2M n J (2tt) 4 \dk 



D(k s 



(4.6) 



tr 



j(-k,P)S(k q )^(k,P) 



N x = ~^r^i 



d 4 p d 4 k 1 



2M n J (2tt) 4 (2tt) 4 2iV 2 



P{- P l)P{-pl) (4.7) 



tr 



i>(-p,P) 



d 



S(q))Mk,P) 



Np = --r^i 



d A p d A k 1 



2M„ " J (2tt) 4 (2tt) 4 2iV 2 

p^P'i-pDPi-pl) - p 2 ,P(-p\)P'(-pl) 



(4.8) 



tr 



i>(-p,P)S(q)Tp(k,P) 



withpi = -(l-3»/)P/2+p + fc/2, p 2 = (l-3r?)P/2- 
p/2-k, for (j = a' = l/2. 



It is clear from Eq. (4.4) that either the exchanged 
quark or the presence of a diquark substructure, or in 
general both, provide non-vanishing contributions to the 
nucleon BS normalization. Maintaining these terms is crit- 
ical for the correct determination of the electromagnetic 
charges of baryons since the nucleon normalization and 
electromagnetic form factors are intimately related by the 
differential Ward identities. 

In the following section, we show that in calculations 
of the electromagnetic form factors of the nucleon, use of 
the usual impulse approximation p3| (which includes only 
contributions arising from the coupling of the photon to 
the quark and diquark, that are related to the terms N q 
and No in the normalization condition), by itself is insuf- 
ficient to guarantee electromagnetic current and charge 
conservation of the nucleon. 

One of the additional contributions required for the 
proper calculation of electromagnetic form factors that 
goes beyond the usual impulse approximation is the cou- 
pling of the photon to the exchanged-quark in the kernel 
from Eq. (|3.5|). We will show that this term provides a 
crucial contribution to the nucleon electromagnetic form 
factors and helps maintain electromagnetic current con- 
servation of the nucleon. It is interesting to note that the 
contribution of this term is important even in the spe- 
cial case that 77 = 1/2 in which the term Nx d oes not 
contribute to the normalization condition of Eq. (4.4). 

We will show in the next section that in the presence 
of a diquark with a non-trivial substructure, additional, 
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direct couplings of the photon to this substructure are 
required to maintain current conservation. Similar contri- 
butions have previously been found important in several 
other contexts |44|,|^|]. These couplings are commonly re- 
ferred to as "seagulls" . 

We conclude this section by summarizing that the re- 
spective con trib utions to the n ucleo n n orm alization condi- 
tion of Eq. (4.4) given in Eqs. (4.5) to ( O ) correspond to 
the usual impulse approximate contributions arising from 
the constituent quark and diquark, plus the contribution 
from the exchange-quark in the BS kernel, and the seagull 
contributions which arise from the quark-substructure in 
the diquark correlations. 



5 The Electromagnetic Current in Impulse 
Approximation and Beyond 



The electromagnetic current operator J™ (a;) in impulse 
approximation is determined by the disconnected contri- 
butions from the electromagnetic couplings of the specta- 
tor quark (Jjf) and the scalar diquark (J^) which in the 
Mandelstam formalism are calculated from the following 
momentum space kernels J43|, 



J^(p,P';k,P) = {2n)H\p-k-f l Q) 

q g r?(r,P' +p, V P + k)D- 1 (r]P-k) , 
j£(p,P';k,P) = (27T) 4 S 4 (p-k + V Q) 

q D p£ (fjP' -p,fjP- k) s-^rjP + k) , 



(5.1) 



(5.2) 



with Q = P' — P and j]+fj = 1. The charge of the spectator 
quark in the nucleon is denoted q q , the charge of the scalar 
diquark is qo, and q q + qo = 1 and for the proton and 
neutron, respectively. The Ward-Takahashi identities for 
the quark and diquark electromagnetic vertices are 



Q^(p + Q, P ) = iS-\p + Q) - iS~\p) , 



Q»r£(p + Q,p) = iD-\p + Q) - iD-^p) 



(5.3) 
(5.4) 



From these one can immediately write down the nucleon 
matrix elements for the divergences of the correspond- 
ing Mandelstam currents between initial and final nucleon 
states with momentum and spin P, s and P', s' respec- 
tively, as 

(P',s'\d^(0)\P,s)= q q J (5.5) 

u(P', s') i>{-{k + fjQ), P')i>{k, P) u{P, s) 

- u(P', s') + fjQ), P')if(k, P) u(P, s)}, 

(P',s'\d^(0)\P,s) = to / (5-6) 

u(P', s') $(-(k - rjQ), P')ip(k, P) u(P, s) 

~ u(P', s') i>{-{k - V Q),P')$(k, P) u(P, s)} . 



Here, we insert the BSEs for the amplitudes ip and ip which 
can be written in the compact form, 



d*k 
d 4 k 



K(p,k,P)ip(p,P) , 



ij(-k,P)K(p,k,P). 



(5.7) 
(5.8) 



After shifting the integratio n mo ment um by p + fjQ — > p 
in the second terms of Eqs. fl5.5|) and (5.6), one obtains, 



{p'^jm^) = ^j^^{ 

u(P', s') i>(-p, P')K{p, k + fjQ, P')ip(k, P) u{P, s) 

- u(P', s') 4>{-p, P')K{p - f,Q, k, P)<P(k, P) u{P, s)} , 

(P',s'\dM0)\P,s) = q D J (5.10) 

u(P', s') j(-p, P')K(p, k - V Q, P'Mk, P) u(P, s) 

- u(P', s') ip(-p, P')K{p + V Q, k, P)^{k, P) u(P, s)\ . 



(5.9) 



Examination of the terms in Eqs. (5J3) and ( 5.10 ) reveals 
that their sum gives rise to a conserved current only if 
K(p, k, P) = K(p — k). That is, if the nucleon BS kernel 
is independent of the total nucleon bound-state momen- 
tum P and, in addition, it only depends on the difference 
of the relative momenta p — k. These criteria are satis- 
fied in studies of meson bound states within the ladder 



approximation, see for example Rcf. 34 1. However, we ob- 
serve that even in absence of an explicit dependence on 
the total nucleon momentum P, the exchange kernel of 
the BSE, as obtained from the nucleon Faddeev equation, 
necessarily depends on the sum of the relative momenta 
p + k and not their difference. It follows that even with 
approximating the exchange kernel by a P-independent 
one, which corresponds to neglecting diquark substruc- 
ture together with using 77 = 1/2 as in Refs. |3(||3^J, the 
electromagnetic current of the nucleon is not conserved in 
the impulse approximation and it is already necessary to 
include an additional coupling of the photon to the quark- 
exchange kernel. 

It is interesting to note that similar photon-kernel cou- 
plings are required in other systems as well. For example, 
while the nucleon-nucleon scattering kernels due to meson 
exchanges depend only on the difference of the relative 
momenta (i.e., K = K(p — k)), the isospin dependence 
of a charged meson that is exchanged between the two 
nucleons requires one to introduce additional photon cou- 
plings to the exchanged meson to maintain current conser- 
vation, in this case the charged pion. Such contributions 
play an important role in determining the electromagnetic 
form factors of the deuteron ^(| . For the importance of 
meson-exchange currents in few-nucleon systems, see also 
the recent review, Ref. |47j]. 

The coupling of the exchange qua rk ( with electromag- 
netic charge qx) in the kernel of Eq. (3J3) gives rise to the 
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additional contribution Jx to the nucleon current: 

J x (p,P';k,P) = -Qx\ (5.11) 

X(pi ,fjP~k) S T (q) rf (q' , q) S T (q') f (p' 2 , fjP' - p) , 
with q = f\P — r\P' — p — k , q' = q + Q , 
p x = a{r,P' +p) -aq, p' 2 = -a\r,P + k) + a'q' , 
and <7 + <7 = <7+<t=1. 

From the War d- Ta kahashi identity for the quark-photon 
vertex in Eq. (5.3), one finds that the divergence of the 
exchanged-quark contribution to the nucleon electromag- 
netic current is 



s „J x (p,P';k,P) = 



(5.12) 



Mx ~ ( xiPuVP ~ k ) S T (q) X (p'z, f/P' - P) 



XipufjP - k) S T (q') fip'^fjP' -p)} 



To provide a comparison with the quark and diquark elec- 
tromagnetic curr ents given above for the quark-exchange 
kernel of Eq. ( |3.5| ) writte n usi ng the same momentum con- 
ventions as used in Eq. ( 5.11 ), we rewrite the diverg ences 
of the Mandelstam currents given in Eqs. ( |5.9| ) and (5.10) 
as 



Q^{p,P';k,P) = 



(5.13) 



\ (x(pi,flP-k)S T (q) X T (P2-Q,vP' -P) 



X(pi -Q,fjP- k) S 1 (q') X (p' 2 , VP' - P) 



QM Pl P';k,P) = 
1 



(5.14) 



"to 2 \ x{pi-vQ,fiP-k + Q)S T (q') X (p^vP'-p) 

- X(P1,VP - k) S T (q) f{p' 2 - fr'Q, rjP'-p-Q)) . 

Since q q — qo + qx = 0, one thus finds that, in the case 
of a point-like diquark BS amplitude, (i.e., neglecting any 
momentum dependence in the diquark BS amplitudes x 
and x)i the sum of the three currents given above now 
yields a conserved electromagnetic current for the nucleon; 
that is, 



■I' \ J q 



■'d 



1*' 

• J x 







(5.15) 



For this cancellation all three contributions are crucial. In 
particular, the contributions from the photon coupling to 
the exchanged quark J x , as well as the impulse approxi- 
mate contributions + must all be included. For the 
quark-diquark model of baryons, these three contributions 
to the current correspond to those given in Ref. Jl8| for 
the general structure of the one-particle contributions to 
the current of a 3-particle Faddeev bound state (when 
the interactions are due to a separable 2-particle scatter- 
ing kernel)^] In the NJL model the one-particle contribu- 

4 Here, one-particle refers to a contribution that arises from a 
one-particle irreducible 3-point vertex for the photon coupling. 



tions 
state 



,deld the complete current of the 3-particle bound 



However, if the substructure of the diquark BS am- 
plitudes is taken into account and the diquark BS am- 
plitude is dependent on any momentum, the cancellation 
of the long itudi nal pieces in the one-particle contribu- 
tions, Eqs. ( 5.12| ) to (5.14), is destroyed. Additional pho- 
ton couplings which are not of the one-particle type be- 
come necessary. The violations to current conservation 
from the one-particle contributions can be displayed in 
the present framework in a way which will become useful 
in foll owin g s ectio ns. Re arran ging the six terms from the 
Eqs. ( 5.12|) , ( 5.13| ) and ( |5.14 ) in such a way as to factor 
out two terms, 

Si{p,P';k,P) := -iq q x( Pl -Q,f)P-k) + (5.16) 
iQD x(pi - o-Q, f]P - k + Q) - iq x x(Pi,vP - k) , 



S 2 {p,P';k,P) := -iq q x(P2- 
ilDZip'?, - o-'Q,f]P' -p 
one obtains 



Q 1 f j P / -p)+ (5.17) 
- Q) - iqx^(p' 2 ,vP' ~P) , 



Ql 



J!, 1 



J x ) = 



(5.18) 



S 1 (p,P';k,P)S T (q')x (p' 2 ,f)P'-p) 



~x(PuVP-k)S T (q)S^(p,P';k,P) 

It will be demonstrated in the subsection below that these 
contributions are exactly canceled by the so-called "seag- 
ull" contributions which arise from one-particle irreducible 
4-point couplings of the two quarks, the diquark and pho- 
ton. Such terms must be included whenever the substruc- 
ture of a diquark bound state (in form of a momentum- 
dependent diquark BS amplitude) is included in the de- 
scription of the nucleon. Analogous seagull contributions 
were previously found necessary in 7-meson-baryon-baryon 
couplings to satisfy the corresponding Ward- Takahashi 
identities f44,||j|. 



Note that terms analogous to the explicit one-particle 
contributions to the bound state currents presented in 
this section can also be obtained by employing a "gen- 
eralized" impulse approximation for the 3-particle Fad- 
deev amplitudes. This procedure was recently adopted in 
an exploratory study of the electromagnetic nucleon form 
factors Q. In this study five distinct (one-particle) con- 
tributions to the form factors arose which were calculated 
using parameterizations of a simplified nucleon Faddeev 
amplitude. From the separable-kernel Faddeev equation 
one readily verifies, however, that only three of these five 
contributions are independent. These three have the ex- 
act same topology as the one-particle contributions pre- 
sented above. Starting from the generalized impulse ap- 
proximation, the relative weights of these contributions 
differ, however, from those needed for current conserva- 
tion. The latter can be systematically derived from a gaug- 
ing technique |50 . The discrepancy in the weights is due 



to an overcounting of the (generalized) impulse approxi- 
mation which can therefore not lead to a conserved cur- 
rent plj . This problem is independent of the necessity for 
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the additional seagull contributions which persists when 
non-pointlikc diquarks are used. We will now address these 
contributions. 



5.1 Ward Identities and Seagulls 

The Ward-Takahashi identity for the quark-photon vertex, 
Eq. ( |5.3| ), follows from the equal- time commutation rela- 
tion for the electromagnetic quark-current operator jy,(x) 
with the quark field (with charge q q ), 

[f(x),q(y)]S(x Q -y ) = -q q q(x)5 4 (x ~ y) , 

\j°(x),q(y)}S(x -y ) = q q q(x) S 4 (x - y) . (5.19) 

Formal problems with equal-time commutation relations 
for interacting fields can be avoided by replacing the canon- 
ical formalism with a Lagrangian formulation based on rel- 
ativistic causality rather than to single out a sharp time- 
like surface |^2| . However, as an operational device for the 
derivation of War d iden tities, the equal-time commutation 
relations of Eqs. ( |5.19 ) will nevertheless give the correct 
result. 

Consider the 5-point Green function that describes the 
photon coupling to four quarks G^ 7 g S . Using the notation 
that q qa , q q p, q ql , and q q g denote the charges of the quark 
fields with Dirac indices denoted by a, 0, 7 and 6, respec- 
tively, the Ward identity is given by 

dl b {T{q 1 {x i )q a {x 1 )qp{x 2 )q s {x i )f{z))) = (5.20) 

- ( q qa 5 4 {xi - z) + q ql S 4 (x 3 - z) - q q/3 5 4 (x 2 - z) 

~ q qS 6 4 (x 4 - z)) {T{q 1 {xz)q a (x 1 )qp{x2)qs{xi))) ■ 

The 4-point function on the right-h and side has the di- 
quark pole contribution given in Eq. (2.2) and depicted in 
Fig. B. Th e Fourier transformation of the left-hand side of 
Eq. ( p>.20D allows one to define, 

(5.21) 



d 4 xi d 4 x 2 d 4 x 3 d 4 x 4 e lpaX1 e' p ' ,X2 e lp - X3 e lpsXi 

(T(q~ t (x 3 )q a (x 1 )q l3 (x2)qs(x4)j tJ '(0))) ■ 

Here, p = ap 1 — ap a , k = cr'pp — cr'ps and P' = P + Q 
as before. It is straight forward to verify the following 
Ward identity for this 5-point Green function from the 
pole co ntri bution to the 4-quark Green function, given 
in Eq. (2.2), which determines the dominant contribution 
when the diquark momenta are close to the diquark pole 
at P 2 = P' 2 = ml. 



%3 



q q s 



a pi 


— m 2 , 
i 


+ ie 


p2_ 


- m 2 -+ 
i 


. A 

it 


p/2 


i 


He 



p, 



(5.22) 

X ia (p + ctQ,P) X(3s{k, P) 
X~fa(p- crQ,P) Xf3s(k,P) 



X^a(p,P')Xi3s(k + <j'Q,P') 



To explicitly demonstrate that this does indeed give the 
additional contributions necessary to current conservation 
of the BSE solution for the nucleon, one needs to con- 
sider the irreducible 4-point coupling of the photon to the 
quarks and the diquark derived from the following defini- 
tion: 



(5(p 7 ) M"fa,p a ,P d ) S T (p a )) ia D(P d ) 
d 4 k 



(5.23) 

GJL Pd + Q; k, P d ) xsp(k, P d ) , 



(2tt) 4 Q7 ' /3 ' 5 

with p a = -p + <r(P d + Q), p y =p+ cr(P d + Q), and 
d 4 k 



Z := 



(27T)' 



tr \x(k,P d )x(k,P d ) 



(5.24) 



The Ward identity for the 5-point Green function then 
entails, 

iQ^{p llPa ,P d ) = (5.25) 

Qqa X(P + vQ,Pd) S T (p a - Q) S~ 1T (p a ) 

+ q g p S^ 1 {p 1 ) S(p 7 - Q) x(p - crQ, P d ) 



-A*(Q 2 )x(p,P d + Q) 



Pj - ml 



(Pd + Q) 2 -ml+ie' 



with 



A$(Q 2 



Z 



(5.26) 



Q = p 7 + p a - P d , p = ap 7 - apa , 
(<r + a = 1) and A<t,(Q 2 ) is defined by 
-1 f d 4 k r 

J (27T)4\ 

q qP tr[S T (-k + a'P d + Q) x 

X(k - a'Q, P d + Q) S(k + <j'P d ) x(k, P d ) } 
+ q qS tr[S T (-k + cr'P d ) X 

l{k + a'Q, P d + Q) S(k + a'P d + Q) x{k, P d )] } . 

In the limit Q — > 0, this is normalized in such a way as 
to yield the charge of the scalar diquark; that is, A$(0) = 

q q @ + q q s = q<s>- 

In a more detailed and complete calculation, the cou- 
pling of the diquark to the photon would itself have to 
be done within a Mandels tam formalism. To achieve this, 
the Ward identity of Eq. (|5.3| ) would be used in the Man- 
delstam formalism to construct the Ward identity for the 
quark substructure of the diq uark , thereby replacing the 
naive Ward identity of Eq. ( |5.4[ ) with a more accurate 
identity which accurately depicts the quark substructure 
of the diquark. This added complication can be worked 
out in a straight-forward manner by introducing a few 
additional technical details. However, the basic principle 
of such couplings, as derived from Ward identities, can be 
seen from the following simplifying assumption, which will 
be used in the following sections. Assume that A^> (Q 2 ) is 
independent of the photon momentum, such that A$(Q 2 ) = 
A$(0) — q&. This assumption is sufficient in order to ob- 
tain the correct charges for the nucleon bound state. From 
this starting point, the electromagnetic diquark form fac- 
tor can be easily included in a minor extension of the 
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framework and follows simply from the inclusion of a de- 
pendence on the photon momentum Q 2 of A${Q 2 ). 

By including the effect of only the charge of the di- 
quark, that is setting A^Q 2 ) = q® for all photon mo- 
menta Q, the divergence of the amplitude is written 
as 

Q^AP( Pa ,p fj ,P d ) = (5.27) 

(FM^'btoPfrPa) + Q»MF(p a , Pfl ,P d ) , 

where M legs contains the couplings of the photon to the 
amputated quark and diquark legs according to their re- 
spective Ward identities, 

iQ^M l ^ s (p a ,p ,P d ) = (5.28) 
q qa (S- 1 ( Pa )-S- 1 (p a -Q)) x 

S(p a -Q)x(p-*Q,Pd) 
+ q ql 3x{p + (tQ, P d )s T {pp -Q)x 

(s-^m-s-' t ( P0 -q)) 

- q$x(p, Pd + Q)D(P d + Q) 

(D-^P^-D-^Pa + Q)) , 

where the term M sg describes the one-particle irreducible 
seagull couplings and its divergence is given by 

iQ»M?(p a ,P(i,Pd) = (5.29) 
q qa x(p - Pd) + q q/ 3 x(p + (?Q, Pd) 
- q$x{p, Pd + Q) ■ 

These seagull couplings are exactly what is needed to 
arrive at a conserved electromagnetic current for the nu- 
cleon. Upon substitution of the charges of the spectator 
and exchanged quark and the scalar diquark, q qa = q q , 
qqp — qx and q$ = qn, respectively, one finds 

S^P'-k.P) = Qni s ^(r,P' +p,q + Q lV P-k), (5.30) 

with Q = P' — P. A solution to this Ward identity, with 
transverse terms added so as to keep the limit Q — ► 
regular, which follows from a standard construction, c./., 
Refs. @f§j|, is provided by 

iM^( Pa , Pf3 ,P d ) = (5.31) 
(2p a - Q) M 

'pi - (p« - QY 
q p}-(p0-Q) 2 



'''!~T~ 77, 777)2 1 X(P a ~ Q'PPi P d) ~ X.{Pa,Pt3, Pd) 

' X(Pa,Pl3-Q,Pd)-x(Pa,Pl3,P d ) 



qD- 



(2P d + Q\ 



X{Pa,Pf3,Pd + Q) ~ x(Pa,P/3,Pd) 



'{P d + Q) 2 -P 2 d 

with Q = p a +pp — Pd- Analogously, for the other seagull 
one finds 



iM°9( Pa , Pp ,P d ) = 
(2p a - Q)jj_ 

Q 

(2 Pg -Q) M 

qx pl-{p0-Q) 2 



(5.32) 



q<r 



{p a - QY 



X(Pa ~ Q,P0,Pd) - X(Pa,Pp,Pd) 
^(Pa,P/3 ~ Q,Pd) ~ X(Pa,P/3,Pd) 



qD 



' X(P*,Pf3,Pd- Q) - x(Pa,P(3,Pd) 



P 2 d -{P d -QY 



The amplitudes x(Pa,P/3, Pd) herein need to be constructed 
from the BS amplitude of the scalar diquark by removing 
the overall momentum conserving constraint p a +Pfj = Pd- 
With these seagull couplings, a conserved electromagnetic 
current operator is obtained by including the seagull con- 
tribution 



J^(p,P;k,P) 



(5.33) 



M s J>{ n P + p, q', fjP - k) S r (q') X (pa, f,P> - p) 



- I X .. 

2 V " 

-x(pi,f)P- k) S T {q) M° gT (-(r,P + k), -q, f]P' - p) 



The total conserved electromagnetic current of the nu- 
cleon is therefore given by J£ m := J£ + + + J£ g - 

Note that this explicit construction of the conserved 
current complies with the general gauging formalism pre- 
sented in Ref. Eq]. In the reduction of Faddeev equa- 
tions with separable 2-particle interactions, the seagull 
couplings arise from 2-particle contributions to the bound 
state current. These contributions describe the irreducible 
coupling of the photon to the 2-particle scattering kernel. 

Technically, the additional contributions to the elec- 
tromagnetic nucleon current arising from exchanged quark 
in the nucleon BSE kernel , as well as the seagull term 
Jjfg, involve two 4-dimensional loop integrations to calcu- 
late the electromagnetic form factors from the nucleon BS 
amplitudes. As demonstrated above, this considerable ex- 
tension to the Mandelstam formalism (involving only sin- 
gle loop integrations) is absolutely necessary to correctly 
include the non-trivial substructure of the diquark corre- 
lations and maintain current conservation of the nucleon. 

While seagull contributions are not necessary in an 
approximation that employs point-like diquarks, as is the 
case in Ref. [ p6[ , beyond-impulse contributions, such as 
the coupling of the exchangcd-quark to the photon are 
necessary! In the study of Ref. Q , it was observed that 
neglecting this contribution produced negligible violations 
to the charges of proton and neutron and so it was dis- 
missed as unimportant. However, in the case of the present 
study this contribution is significant. The reason for this 
is the larger value of the coupling strength g s (obtained 
from the diquark normalization condition g s = \/N 2 ) 
used herein. For a point-like diquark, the coupling need 
not be as large. Furthermore, in contrast to the impulse- 
approximate Mandelstam currents, which are independent 
of g s , the contribution to the electromagnetic current due 
t o th e exchanged quark is proportional to g 2 (see Eq. 



( 5.11 )). Hence, use of a smaller coupling strength gj in 



the BSE, reduces the importance of going beyond the im- 
pulse approximation. 

As an example of the importance of the contributions 
beyond the impulse terms of J q and Jp> , we consider the re- 
sults for the proton and neutron electromagnetic charges 
using the amplitudes plotted in Fig. |3.2| from the Man- 
delstam current Jjf + alone. This leads to charges 
Qp = 0.85 for the proton and Qm = 0.15 for the neu- 
tron. A theorem constrains the charges of the proton and 
neutrons to obey Qp + Qn — 1- The theorem relies on 
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using 7/ = 1/3 and is deriv ed from the nucleon normaliza- 
tion condition in Eq. (4.4). However, the way it is realized 
here is not very satisfying. In Sec. |6|, we discuss in detail 
the relevance of the various contributions to the electro- 
magnetic form factors of the proton and neutron due to 
exchanged-quark-photon coupling and seagull couplings. 

Having proven that the present framework conserves 
the electromagnetic current, one might think that the re- 
sults for the proton and neutron charges, Qp = 1 and 
Qn = 0, must follow trivially. However, the verification 
that this framework provides the correct charges for the 
proton and neutron is not entirely trivial as is demon- 
strated in the next section. 

For finite momentum transfer Q = P'—P > 0, the com- 
plete Mandelstam couplings for the diquark will still re- 
quire modifications. If the photon is coupled to the elemen- 
tary carriers of charge only, that is, to the quarks within 
the diquark (with cjusxk charges Qqa , q q p), the photon- 
diquark vertex will itself be of the form, 



(5.34) 



d^k 

(2tt)« 



tr 



S T (-k + aP d + Q) x 



X(k - *Q, P d + Q)S(k + aP d ) x(k, P d ) x 
S T (-k + aP d )r f " T (-k + aP d + Q, -k + aP d ) 

-14 I 



z q ^Jwr tT 



S 1 (-k + aP d ) x 



X (k + &'Q, P d + Q)S(k + aP d + Q)x 
r»{k + aP d + Q,k + aP d ) S(k + aP d ) x{k, P d ) 



This gives the correct diquark charge in the limit the pho- 
ton momentum Q — > 0. In this limit, as far as the electric 
form factors of the nucleons are concerned this detail is 
irrelevant since they are constrained to be proportional to 
the charge of the nucleon. On the other hand, the anoma- 
lous magnetic moments of the nucleons may receive ad- 
ditional contributions from the quark substructure of the 
diquark. 



6 Electromagnetic Form Factors of the 
Nucleon 



The matrix elements of the nucleon current can be para- 
metrised as 



(P',s'\J^(0)\P,s) 
u(P',s') 



(6.1) 



7^ 1 + |gV^ 



u(P, a) 



Here, T\ and Ti are the Dirac charge and the Pauli anoma- 
lous magnetic form factors, respectively |5j| . Q M = P' — P 
is the spacelike momentum of the virtual photon probing 



the nucleon (— Q > 0). Using the Gordon decomposition 



u{P',s') 



2M 



-u(P,s) 



u(P',s') 



r 



J BF 



u(P, s) , 



(6.2) 



(P> 



with the definition of the Breit momentum P BF 
P)/2, the current can be rewritten as 

(P',s'\J? m (0)\P,s)= (6.3) 



u(P',s') 



r (^i + kt 2 ) 



dp 

r B 



u{P lS ) . 



It is convenient in the following to introduce (matrix val- 
ued) matrix elements by initial and final spin-summations, 

(P'|>|P) := (P\s'\J»\P,s) £u(PV)u(P,«) , (6.4) 



to remove the nucleon spinors. The frequently used Sachs 
electric and magnetic form factors Ge and Gm, are intro- 
duced via 

r>2 



° E ~ Tl + 4M2 KJC " 2 ' 
Gm = T\ + kTi ■ 

These can be extracted from Eq 

(p'\JL(o)\p) = 

A+(P') 



(6.5) 
(6.6) 



^Gm + M -0- {Ge-Gm 

BF 



by taking traces of ( 
M 



Gi 



2p2 ' ^ (^/i ^BF 1 



(6.7) 
A+(P), 

P'\JL(0)\P) as follows: 

(6.8) 



BF 

r2 



M* ( M 
Gm = i tr {Jp-f - pr tr < J P P & 



(6.9) 



We calculate the current matrix elements using Mandel- 
stam's approach with the current operators defined in the 
previous sections, such that 

(p'\JL(o)\p) = 

d 4 p d 4 k 



(6.10) 

(2tt)4 (2n) 4 ^ P ' ] J ^ P ' ; fc ' P) ^ P) 

The current operator consists of the four parts which 
describe the coupling of the photon to quark or diquark, to 
the exchanged-quark and the seagull contributions which 
arise from the coupling of the photon to the diquark BS 
amplitudes. These are determined by the following kernels, 

q q r»(p q ,k q ) D-\k s ) (2n) 4 5 4 (p -k-rjQ), (6.11) 

q D r£(p s ,k s ) S- 1 (k q )(2n) 4 S' i (p~k + f 1 Q) , (6.12) 

-Qxl (6.13) 



J? 



X( Pl ,k s ) S T (q) rftf, q) S T {q') X (p' 2 ,Ps) 
l -{M^ g {p q ,q',k s )S T {q l )f{ J) ' 2l p s ) 



(6.14) 



x( Pl ,k s )S T {q)M» g T (-k q ,-q,p s )) 
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-fc S +Q , Pg— k q 



p s —k s , p q — k q +Q. 

Fig. 6.1. Impulse approximation diagrams. 








k 9' 


r 




T 




v 



q=(f/-^i)P BF -p-k-Q/2, 
q'=(fj-ri)P BF - P -k + Q/2, 
Pl=(^l)Q/4-(l-3»))P BF /2+P+fc/2, 
p'l = (»7-l)Q/4-(l-377)PBF/2+P+fc/2, 
P2=(i?-l)Q/4+(l-3^)P BF /2-p/2-fc, 
p^ = (TH-l)Q/4+(l-3T,)P BF /2-p/2-fc. 




The abbreviations for the various momenta are summa- 
rized in the following table: 





incoming 


outgoing 


quark: 




=r]P+k 


Pg = 


=,,p' 


+p 


diquark: 




= ?jP-fc 


Ps = 


=f,p' 


-p 


exchange quark: 


9= 


fjP—riP'—p—k 


q'= 


=ijp' 


— 7]P— p — k 


relative momenta 


(- 


=<r'=l/2) 








within diquark: 


Pi= 


= jf (Pq-l) 


P2 = 




-kg+q') 


seagull quark-pair: 


p'i- 


= j (Pq-l') 


P2 = 


=§(- 


-k q +q) 



The contributions to the form fac tors in the impulse 
approximation are depicted in Fig. 3.1 , while the exchange- 



quark and seagull contributions are shown in Fig. 6.2. 

We use a = a' = 1/2 in the diquark amplitudes. As 
discussed in Sec. ^, this implies that the relative momenta 
Pi within the diquarks are exchange symmetric and that 
our parameterizations of the diquark BS amplitudes are 
independent of the mass of the diquark. We thus set, 



N s 



(6.16) 



Php'i) 



This also simplifies the seagull terms, as the seagull cou- 
plings to the diquark legs do not contribute to the seagulls 
in this c ase. T he am plitu des in the brackets of the last lines 
of Eqs. (|573l| ) and (|||) cancel. 

The construction of such vertex functions from the 
Ward-Takahashi identities is not unique. In particular, the 
forms for the irreducible seagull couplings Mff q and Mjf q 
given in Eqs. ( 5.3lj) and ( 5.32 ), respectively, are designed 
for amplitudes \ an d X which are functions of the scalars 
Pai P% an d' m general, Pj. For Q — > 0, the possibility that 
the denominators in each of the three terms may vanish 
entails that the prefactors that arise from expanding the 
amplitudes in brackets must also vanish. 



x((p*-Q) 2 ,pD 



x(p 2 a ,p}) 



(6.17) 



2 (P»Q) QjpXiPaiPp) 



Our present assumption on the dominant momentum de- 
pendence of these amplitudes is slightly d iffere nt though. 
For amplitudes x = x{{Pa — Pp) 2 /4), see ( 6.16 ), the pref- 
actors of the seagulls in the form given in Eqs. ( 5.31 ) 



Fig. 6.2. Exchange quark and seagull diagrams. 



and ( |5.32j ), corresponding to factors oc \/{p a Q) and oc 
^/{ppQ) respectively, are not canceled in an analogous 
way. To cure this, we replace the quark momenta p a and 
pp in these prefactors by (plus/minus) the relative mo- 
mentum, ±(p Q —pp)/2. This yields, 



MS, 



l M P L - Q/2) 2 ) - X(P?)) (6.18) 



Ap[Q - Q2 

(M + QY 



-qx 



4p[Q + Q 2 

(4p 2 - QY 



(x(ip[ + Qm 2 )-x(p' 2 )) 



4 P2 g - q 2 

(4p 2 + QY 



+qx 



4 P2 Q + Q 2 



(x((p 2 -Q/2) 2 )-x(^)) (6.19) 

mp2+Q/2) 2 )-w2)) 



Since Ward-Takahashi identities do not completely con- 
strain the form of the vertices, such modification of this 
type is within th e freedom allowed by this ambiguity. The 
forms ( 6.18 , 6.19 ) solv e the corresponding Ward identities 
at finite Q, see Eq. ( 5.29j ) (with q q + qx = to)- In ad- 
dition, the smoothness of the limit Q — > for the given 
model assumptions is ensured. In this limit, 



{q q - qx) Pi x'ipl) , 
{q q - qx) P2 x{p\) ■ 



(6.20) 



We emphasize that this limit is unambiguous. It must co- 
incide with the form required by the differential form of 
the Ward identity for the seagull couplings. As such it re- 
stricts contributions that are both longitudinal and trans- 
verse to the photon four- momentum Q^. The form above 
follows necessarily for the model diquark amplitudes em- 
ployed in the present study, and this form provides the cru- 
cial condition on the seagull couplings that ensure charge 
conservation for the nucleon bound state. 

The nucleon charges are obtained by calculating 



d A p d 4 k 



(6.21) 



4 (2tt) 4 

tr$(-p, P) pMp, P; k, P) V>(fc, P)} 
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The various contributions to the electromagnetic current 
for Q — ► {i.e., P' = P) are given by 



D 



J sg 



u hl j JLs-\k q fj D-\k s ) {2^5\p - k), (6.22) 
iqn (JkJ D ~ 1{ks) ) S ~ 1{kq) ^ 5 ^P- k ^ ( 6 - 23 ) 

iqx 2^ p (-*?) p H# (a^) ' (6 - 24) 



1 



27V? 



S(<z) 



(6.25) 



{pi,p\-p\)p{-pI) - p^P{-p\)P'{-pD) . 

Comparing this to the normalization integrals given in 
Sec. [|, one finds 

G E (0) = (6.26) 
q q N q + q D N D + q x N x - (q g - qx) N P . 

Using q q — 2/3, qo = 1/3, qx = — 1/3 for the proton, and 
q q = —1/ 3, gp = 1/3, qx = 2/3 for the neutron, together 
with Eq. ( |4.4| ), one therefore has, 

1 = r]N g + (1 - ?7)iV D + (1 - 2 V )N X + (1 - 3tj)JV P , 
Qp = \N q + ^N D - ±N X - N P , (6.27) 



Qn = -3^9^ 3 



1 2 

N D + ^N x 



N, 



p ■ 



However, these three equations are not independent. Rewrit- 
ing the normalization condition for the nucleon BS ampli- 
tudes, we find that 

1 = l Nq+ l N ° ~ l Nx ~ Np (6 - 28) 

+ {V ~ 5) (N q -N D - 2N x - 3N P ) , 
which entails that 

1 = Q P + (2- 3 V )Q N . (6.29) 

To verify that we do in fact obtain the correct charges of 
the proton and neutron, it suffices to show that N q — Np, = 
2Nx + 3Np; that is, it suffices to show that the neutron is 
neutral, Qn = 0. The proof of this is straightforward and 
is given in Appendix |c[ 



6.1 Numerical Computation 

The numerical computation of the form factors is done in 
the Breit frame, where 



Q M - (0,Q) , 
P" = (wq,-Q/2) 
P'" = (w ,Q/2) , 
P& = (« Q ,0) , 



(6.30) 



The transformation of these variables to 4-dimensional 
Euclidean polar coordinates follows the sam e prescriptions 
as those employed in Sec. || (see Eqs. ( 3.19] )), namely 

(6.31) 



{-p 2 , -k 2 , -Q 2 } , Pi 

k |Q| zq , 

iujq k z BF , 



{p\ k 2 , Q 2 } 

pQ -p\Q\vq , kQ - 

pP BF -> iuj Q py BF , kP B 
pP' = P P BF +pQ/2 -> 

iuj Q py BF - p\Q\vq/2 =: iMpy, 
kP = kP BF - kQ/2 -> 

iwq k z BF + k\Q \ zq/2 =: i M k z . 

In the presence of two independent external momenta, Q 
and P BF , we are left with 5 independent angular variables. 
Together with the absolute values of the integration mo- 
menta p and k the exchange-quark and seagull contribu- 
tions to the form factors at finite momentum transfer Q 
require performing 7-dimcnsional integrations. These are 
computed numerically using Monte Carlo integrations. 

For the impulse approximation diagrams, the num- 
ber of necessary integrations collapses to three due t o the 
mome ntum-conserving: delta functions in Eqs. ( 6.11 ) and 
( 3.12 ). One of the integrations is over the absolute value 
of the loop momentum k and two are the angular integra- 
tions over zb f and zq , the cosines of the angles between 
k and Pbf and k and Q, respectively. 

The BS amplitudes for the nucleon bound states are 
given in terms of t he tw o scal ar fu nctions Si (p, P) and 
S 2 (p,P), c.f Eqs. (|3T6|) and (§T7|) in Sec. |, which, we 
recall Eqs. (3.20) and (3.21), are expanded in terms of 
Chcbyshev polynomials T n to account for their depen- 
dence on the azimuthal Euclidean variable, 



N-l 



S{p,y) ~ J2(-i) n S n (p)T n (y) 



(6.32) 



While the argument y of the Chebyshev polynomials, the 
cosine between relative and total momentum, is in [—1, 1] 
in the rest-frame of the nucleon, this cannot be simultane- 
ously true for the corresponding arguments in the initial 
and final nucleon bound-state amplitudes at finite (space- 
like) momentum transfer Q. In the Breit frame, these ar- 
guments are, 



M 
M Vbf 



— 1 — 
2 



1 IQ 



-r-r zq and 



M 



.1 IQI 

l 2M- yQ 



(6.33) 



with LOQ = y/M 2 + Q 2 /4. 



for the initial and final nucleon BS amplitudes respectively 
(with the angular variables zq, yQ and z BF , y BF all in 
[—1,1]). In order to use the nucleon amplitudes computed 
from the BSE in the rest frame, analytical continuation 
into a complex domain is necessary. This can be justified 
for the bound-state BS wave functions ip (with legs at- 
tached). These can be expressed as vacuum expectation 
values of local and almost local operators and we can re- 
sort to the domain of holomorphy of such expectation val- 
ues to continue the relative momenta of the bound-state 
BS wave function tp(p, P) into the 4-dimensional complex 
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Euclidean space necessary for the computation of Breit- 
frame matrix elements from rest-frame nucleon wave func- 
tions. The necessary analyticity properties are manifest 
in the expansion in terms of Chebyshev polynomials with 
complex arguments. 

There are, in general however, singularities associated 
with the constituent propagators attached to the legs of 
the bound state amplitudes, here given by the free particle 
poles on the constituent mass shells. For sufficiently small 
Q 2 these are outside the complex integration domain. For 
larger Q 2 , these singularities enter the integration domain. 
As the general analyticity arguments apply to wave func- 
tions ip rather than the truncated BS amplitudes ip with 2- 
component structure R(p, y), it is advantageous to expand 
these untruncated BS wave functions directly in terms of 
Chebyshev polynomials (introducing moments R n (p)), 



R(p,y) 



N-l 

E 

Tt = 



(-<) n Rn(p)T„(y) 



(6.34) 



and employ the analyticity of Chebyshev polynomials for 
the BS wave function R. This can be written in terms of 
the two Lorentz-invariant functions R\ and R2 by 

V(p, P) = D( Ps )S( Pq )ij(p, P) = (6.35) 
D{ Ps )S{p q )(s 1 {p,P)A+{P) + S 2 {p,P)~{p,P)A + {P)) 

=: Ri(p, P) A+(P) + R 2 (p, P) S(p, P) A+(P) . 

The price one must pay, however, is a considerably slower 
suppression of t he hi gher Chebyshev moments in the ex- 
pansions in Eq. (3.34) for the BS wave functions compared 
to the much faster suppression observed for the truncated 
amplitudes ip. For example, the fourth moments of the 
truncated nucleon amplitudes ip are shown in Fig. 3.2 of 
Sec. § Their magnitudes are less than 2 orders of magni- 
tude smaller than the leading Chebyshev moments. One 
must include up to 8 Chebyshev moments in the expan- 
sion for the untruncated BS wave functions in order to 
achieve a comparable reduction. 

If the truncated BS amplitudes ip are used in the ex- 
pansion, one must account for the singularities of the quark 
and diquark legs explicitly when these enter the integra- 
tion domain for some finite values of Q 2 . A naive transfor- 
matio n to the Euclidean metric, such as the one given by 
Eqs. (|6~3l"l ). is insufficient. Rather the proper treatment 
of these singularities is required when they come into the 
integration domain. 

For the impulse-approximate contributions to the form 
factors we are able to take the corresponding residues 
into account explicitly in the integration. Although, this 
is somewhat involved, it is described in Appendix [S| For 
these contributions, one can compare both procedures and 
verify numerically that they yield the same, unique results. 
This is demonstrated in App. 

We have to resort to the BS wave function expansion 
for calculating the exchange quark and seagull diagrams, 
however. The residue structure entailed by the structure 
singularities in the constituent quark and diquark propa- 
gators is too complicated in these cases (involving the 7- 
dimensional integrations). The weaker suppression of the 



higher Chebyshev moments and the numerical demands 
of the multidimensional integrals thus lead to limitations 
on the accuracy of these contributions at large Q 2 by the 
available computer resources. 

The Dirac alge bra necess ary to compute Ge and Gm, 
according to Eqs. ( |6.8| ) and (5.9), can be implemented di- 
rectly into our numerical routines. We use the moments 
S n (p) or R n {p) obtained from the nucleon BSE as de- 
scribed in Sec. [| which are real scalar functions with 
positive arguments. These functions are computed on a 
one-dimensional grid of varying momenta with typically 
n p = 80 points. Then spline interpolations are used to ob- 
tain the values of these functions at intermediate values. 
The scalar functions S(j>, P) an d R(p , P) a re then easily 
reconstructed from Eqs. ( |6.32 ) and (6.34), respectively. 
Then complex arguments, as given in Eqs. ( |6~33|) , appear 
in the Chebyshev polynomials when the electromagnetic 
form factors are calculated. 

In the results shown below, the 3-dimcnsional inte- 
grations of the impulse approximation diagrams are per- 
formed using Gauss-Legendre or Gauss-Chebyshev quadra- 
tures, while the 7-dimensional integrations necessary for 
the calculation of the exchange-quark and seagull contri- 
butions are carried out by means of stochastic Monte- 
Carlo integrations with 1.5 xlO 7 grid points. We find that 
beyond Q 2 = 3 GeV 2 , numerical errors for the stochas- 
tic integrations becomes larger than 1% of the numerical 
result. This we attribute to the continuation of the Cheby- 
shev polynomials T n (z) to complex values of z as described 
above. 

In addition to the aforementioned complications, there 
is another bound on the value of Q 2 above which the ex- 
change and seagull diagrams can not be evaluated. It is 
due to the singularities in the diquark amplitudes x(Pi) 
and in the exchange quark propagator. The rational n-pole 
forms of the diquark amplitude, P n _ P (p) — (7,1/(7,1 +p 2 )) n 
for example yield the following upper bound, 



Q 2 < 4 



47™ 



(1-3tj)* 



(6.36) 



A free constituent propagator for the exchange quark gives 
the additional constraint, 



Q 2 < 4 



(l-2rj) 2 -M 2 



(6.37) 



It turns out, however, that these bounds on Q 2 are in- 
significant for the model parameters employed in the cal- 
culations described herein. 



6.2 Results 



In Fig. 6.3, we show the electric and magnetic Sachs form 
factors of the p roto n and neutron using the parameter sets 
given in Table 3.1 in Sec. ^ which correspond to a fixed 
value for the diquark mass m s — 0.66 GeV. The nucleon 
amplitudes used in these calculations correspond to those 
shown in Fig. 1A of Sec. ||. The charge radii obtained by 
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Proton Electric Form Factor 

m s =0.66 GeV, P(p): n-pole or exponential 
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Fig. 6.3. Nucleon electric and magnetic form-factors for fixed widths of 5*1,0. The label Gauss, shifted 
refers to a calculation with diquark vertex function Pgau = exp(— (p 2 /7gau — I) 2 )- The parameters 
are listed in Tabs. [O] and 3.2 in Sec. |J| 



using th e mo del forms of the di qua rk BS amplitude given 
are given in Table 



in Table 3.1 



.1 



Examination of the charge radii given in Table 6.1 re- 
veals that the width obtained for the nucleon BS ampli- 
tude is closely correlated with the obtained value for the 
charge radius of the proton. The dipole, quadrupole and 
exponential forms for the diquark BS amplitude, all give 
reasonable values for the charge radius of the proton. The 
accepted value of the proton charge radius is r p ~ 0.85 fm. 
However, the width of 5i,o(p) does not determine the be- 
havior of the form factors away from Q 2 = 0. This is 
especially clear when the exponential and Gaussian forms 
are employed for the diquark BS amplitude. In this case, 
the proton electric form factor ceases to even vaguely re- 
semble the phenomenological dipole fit to experimental 



data over most of the range of Q 2 shown in Fig. 6.3. The 
neutron electric form factor is even more sensitive to the 
functional form of the diquark amplitudes. The square of 
charge radius of the neutron depends strongly on the cho- 
sen form of the diquark BS amplitude. For the exponential 
and Gaussian forms, the obtained value of is close to 
zero, and it is positive when the Gaussian form of the di- 
quark amplitud e is u sed with its peak away from zero (i.e., 
7^ 0). In Fig. |6.3| , we show the form factors that result 
from using the shifted Gaussian form for the diquark BS 
amplitude with xo/Tgau = 1. In fact, the shifted-Gaussian 
form also produces a node in the electric form factor of the 
neutron, for which there is no experimental evidence. We 
conclude that to obtain a realistic description of nucleons, 
one must rule out the use of forms for the diquark BS 
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Fig. 6.4. Comparison of the electric form factors of proton and neutron with the experimental data 
of Refs. Q,^5| and Ref. [^i), respectivel y. H ere, the dipole form of the diquark amplitude was used 
corresponding to the n = 2 result of Fig. 6.3. 



amplitude which peak away from the origin. 

In Figure 3.4 we compare the n = 2 results, obtained 
from employing the diquark BS amplitude of dipole form, 
to the experimental data of Refs. (54|,|55| and Ref. |3(| for 
the electric form factors of proton and neutron, respec- 
tively. In Ref. Q , the neutron Ge is extracted from data 
taken on an unpolarised deuteron target and with employ- 
ing various TV-TV potentials. As pointed out in Ref. Q, due 
to possible systematic errors of this procedure, these re- 
sults should not be over-interpreted. They can serve to 
give us a feeling for the qualitative behavior and rough 
size of the electric form factor of the neutron, however. 
With resembling the phenomenological dipole fit for the 





form of diquark 


7> (fm) 


r 2 n (fm 2 ) 




amplitude P 


(±0.02) 


(±0.02) 


fixed Si.o-width : 


71=1 


0.78 


-0.17 




n=2 


0.82 


-0.14 




n=4 


0.84 


-0.12 




EXP 


0.83 


-0.04 




GAU 


0.92 


0.01 




GAU shifted 


1.03 


0.37 


fixed masses: 


71=1 


0.97 


-0.24 




71 = 2 


0.82 


-0.14 




71=4 


0.75 


-0.03 




EXP 


0.73 


-0.01 



Table 6.1. The electric charge radii for proton and neutron 
for parameter sets having either the <Si,o-width or the quark 
mass fixed. Error estimates come from the uncertainty in the 7- 
dimensional integration. The corresponding experimental val- 



ues are about r 3 
for the neutron. 



0.85fm for the proton and r„ 



-0.12fm 2 



proton fairly well, as seen in Fig. 3^2 , it might not be too 
surprising to discover good agreement also with the exper- 
imental results for the proton. However, with the special 
emphasis of our present study being put on charge conser- 
vation, such compelling agreement also for finite photon 
momentum transfer and over the considerable range of Q 2 
(from up to 3.5 GeV 2 ), seems quite encouraging. Also 
the neutron electric form factor compares reasonably good 
with the data, especially considering that we did deliber- 
ately not put much effort in adjusting the free parameters 
in our present model. 

The obtained magnetic moments, which range from 

0. 95 . . . 1.26 nuclear magnetons for the proton and from 
-0.80 . . . -1.13 nuclear magnetons for the neutron, are too 
small. The accepted values for the proton and neutron 
magnetic moments are 2.79 and -1.91 nuclear magnetons, 
respectively. The essential reasons for this discrepancy are 
summarized as follows: 

1. The next important diquark correlations which should 
be included in the present framework are those of axi- 
alvector diquarks. These are necessary for an extension 
of the quark-diquark model to the decuplet baryons 
[jj7| . The contribution of their magnetic moments to 
the anomalous magnetic moments of the nucleons was 
assessed in Ref. 57 1. There, the NJL model was em- 
ployed to calculate the electromagnetic form factors 
of on-shell diquarks, and their influence on the nu- 
cleon magnetic moments was estimated from an ad- 
ditive diquark-quark picture. The conclusion from this 
study was that including the magnetic moment of the 
axialvector diquark alone did not improve the nucleon 
magnetic moments. The additionally possible transi- 
tions between scalar to axialvector diquarks, however, 
were found to raise them substantially. Whether this 
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Neutron Electric Form Factor 

m =0.58 GeV, m d =0.66 GeV, P(p): n-pole or exponential 
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Fig. 6.5. Nucleon electric form factors for fixed quark and diquark masses. In the left plot, the bottom 
set of curves represent the results obtained for the proton electromagnetic form factor Ge{Q 2 ) and 
the top set of curves depict the ratios of these form factors over the dipole fit (l + Q 2 /(0.84GeV) 2 ) 2 . 



finding persists in the fully relativistic treatment, is 
subject to current investigations |^8| . 
Furthermore, since the axialvector diquark enhances 
the binding (i.e., lowers the coupling g s ), it tends to 
lower the quark mass required to produce the same nu- 
cleon bound state mass and as will be discussed below, 
a smaller quark mass would also serve to improve the 
obtained values for the magnetic moments. 
In our present study, the diquark-photon vertex is that 
of a free scalar particle. The contribution of the corre- 
sponding impulse approximation diagra m to the mag- 
netic moments (in the right panel of Fig. |6.l| ) is small, 
below 0.01 nuclear magnetons, but non- vanishing. Re- 
solving the diquark substructure by coupling the pho- 
ton directly to the quarks within the diquark, and tak- 
ing into account its sub-leading Dirac structure, will 



dress this vertex as discussed (at the end of Sec. 5.1) 
and it might increase that contribution. We do not ex- 
pect the gain to be substantial though. 
3. The consistency requirement is that the strength of the 
coupling g s — 1/N S is given by the normalization of 
the diquark N s . At present, this leads to a rather nar- 
row nucleon BS amplitude in combination with some- 
what large values for constituent quark mass. Both of 
these effects tend to suppress the quark contribution to 
the form factors arising from the impulse-approximate 
terms. If we had assumed a quark mass m q of 450 MeV, 
and artificially increased the width of the diquark BS 
amplitude, the quark diagram alone would easily con- 
tribute 1.3 . . . 1.5 nuclear magnetons to the magnetic 
moment of the proton. That is, a small change to the 
dynamics of the diquark BS amplitude or mass of the 
quark can have a significant impact on the magnetic 



moment of the proton. A similar sensitivity of mag- 
netic moments of vector mesons to the scales in the 
quark propagator and bound state vector-meson BS 
amplitude was also observed in Ref . |5^] . 

The electric form factors obtained using the fixed val- 
ues for quark and diquark masses of m q = 0.58 GeV and 
m s — 0.66 GeV, respectively, (i.e., for fixed values of the 
binding energy) are shown in Fig. p.q . The corre spon ding 
charge radii are given in the right half of Table 6.1. The 



differences between the various diquark amplitude param- 
eterizations considered herein do not lead to such dramatic 
differing behaviors of the nucleon form factors in this case. 
Nevertheless, we observed that use of the exponential form 
of the diquark BS amplitude still produces a proton elec- 
tric form factor that falls off too fast when compared to 
the phenomenological dipole fit. It also results in a tiny 
value for the square of the charge radius of the neutron. 

Finally, we compare the relative importance of the 
various contributions to the form factors that arise from 
the impulse-approximate, exchange-quark, and seagull di- 
agrams. We separately plot each of these contributions 
to t he t otal proton and neutron electric form factors in 
Fig. 3.6 for comparison. For the purposes of comparison, 
we used the parameter set for the dipole form of the di- 
quark BS amplitude. (This form of the diquark BS ampli- 
tude is used because of the excellent description it provides 
for the the electric form factor of the proton.) 

Here, it is interesting to note that the seagull couplings 
contribute up to about 2 /5 of the proton electric form fac- 
tor G£)(Q 2 )! It is clear that in the nucleon, these beyond- 
the-impulse contributions are certainly not negligible. As 
in the previous discussion of the magnetic moments, where 
we observed a suppression of the impulse-approximate di- 
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Fig. 6.6. Contribution of the single diagrams to neutron and proton electric form factors. The di quar k 
amplitude is the dipole. The sum of these contributions (SUM) corresponds to the result of Fig. 3.4 



agrams arising from the narrow momentum distribution of 
the nucleon BS amplitude which was a result of employing 
narrow diquark BS amplitude. Again, we find that employ- 
ing a diquark BS amplitude that is wider in momentum 
space leads to a wider nucleon BS amplitude and there- 
fore implies that the impulse approximate contributions 
are more dominant than the other contributions. This is 
the case when point-like diquark BS amplitudes are em- 
ployed, such as in Ref. ]36| . 

For the neutron, we observe that its electric form fac- 
tor arises from a sum of large terms which strongly cancel 
each other to produce a small effect. This cancellation is 
the explanation for the appearance of "wiggles" in the re- 
sults for the neut ron Ge(Q 2 ) at moderate Q 2 shown in 
Figs. 6.3, and 6.5. The wiggles are artifacts of the numer- 
ical procedure employed. 



7 Conclusions 

We have introduced an extension of the covariant quark- 
diquark model of baryon bound states. The framework 
developed herein allows for the inclusion of finite-sized di- 
quark correlations in the description of the nucleon bound 
state in a manner which preserves electromagnetic current 
conservation for the first time. For such a framework to 
maintain current conservation of the nucleon, it is neces- 
sary to include contributions to the electromagnetic cur- 
rent which arise from the couplings of the photon to the 
quark-exchange kernel of the nucleon BSE. These contri- 
butions are derived from the Ward-Takahashi identities 
of QED and include the coupling of the photon to the 
exchanged quark in the kernel and the photon coupling 
directly to the BS amplitude of the diquark (the so-called 
seagull contributions). It was shown analytically that the 



resulting nucleon current is conserved and these additions 
are sufficient to ensure the framework provides the correct 
proton and neutron charges independent of the details of 
the model parameters. 

To explore the utility of this framework under the 
most simple model assumptions, simple constituent-quark 
and constituent-diquark propagators and one-parameter 
model diquark BS amplitudes were employed for the nu- 
merical application of the framework. The BS amplitude 
for the scalar diquark was parametrized by the leading 
Dirac structure and various forms for the momentum de- 
pendence of the amplitude were investigated. It was shown 
that the antisymmetry of the diquark amplitude under 
quark exchange places tight constraints on the form of 
the diquark BS amplitude and that the incorporation of 
these constraints have the effect of removing much of the 
model dependence of the diquark BS amplitude parame- 
ters from the calculation of the nucleon electromagnetic 
form factors. 

Calculations of the electromagnetic form factors away 
from Q 2 = require that the nucleon BS amplitudes and 
wave functions be boosted. In this Euclidean-space for- 
mulation, this amounts to a continuation of amplitudes 
and wave functions into complex plane. Two procedures 
to account for proper handling of poles arising from the 
constituent particles in the nucleon are described and com- 
pared. The feasibility of each procedure is discussed and 
explored in detail. It is shown analytically and explicitly 
in a numerical calculation that the two approaches pro- 
duce the same results for the electromagnetic form fac- 
tors. Thus, demonstrating that the framework properly 
accounts for the non-trivial analytic structures of the con- 
stituent propagators. 

For the particular choice of simple dynamical mod- 
els explored herein, the masses of the quark and diquark 
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propagators along with the width of the model diquark BS 
amplitude are the only free parameters in the framework. 
The latter width is thereby implicitly determined from 
fixing the nucleon mass. It is shown, for the case of the 
dipole (and possibly quadrupole) form of the diquark BS 
amplitude, that these few parameters are sufficient to pro- 
vide an excellent description of the electric form factors for 
both the neutron and proton. Other forms of the diquark 
BS amplitude that were explored, such as the exponential 
and Gaussian forms, may be ruled out on phenomenolog- 
ical grounds as they lead to nucleon electromagnetic form 
factors which are inconsistent with the experimental data. 

It was found that the framework is at present unable to 
reproduce the nucleon magnetic moments. The calculated 
magnetic moments are smaller than those obtained in ex- 
periment by about 50%. Possible explanations for this are 
the effect of an inclusion of axialvector diquark correla- 
tions, the addition of more complex structures in the di- 
quark amplitudes, and a resolution of the quark substruc- 
ture in the diquark-photon coupling. These improvements 
are the subject of work currently in progress. 

In conclusion, we find that the covariant quark-diquark 
model of the nucleon provides a framework that is suffi- 
ciently rich to describe the electromagnetic properties of 
the nucleon. However, to ensure that the framework sat- 
isfies electromagnetic current conservation one must go 
beyond the usual impulse approximation diagrams and in- 
clude contributions that arise from the photon couplings 
to the nucleon BSE kernel. In a numerical application of 
this framework, it was found that these contributions pro- 
vide a significant part of the electromagnetic form factors 
of the nucleon and can not be neglected. 
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A Conventions for Diquark Amplitudes 

The BS wave functions X (p, P) and xiPi P) of the (scalar) 
diquark bound-state are defined by the matrix elements, 

X a p(x,y,~P) := (q a {x)q p (y)\P+) (A.l) 

X al3 (x,y;P) := (P + \q a (x)qp(y)) , 

= ( 7 oX t (2/,a;;P)7o) Q/3 (A.2) 

Note that there is no need for time ordering here in con- 
trast to quark-antiquark bound states. The following nor- 
malization of the states is used, 

(P±\P±) = 2w P (2tt) 3 £ 3 (P' - P) , u? P = P 2 + ml , (A.3) 

and the charge conjugate bound state being |P_) = C\P+). 
The contributio n of the charge conjugate bound state is 
included in Eq. (2.2) for Pq — -up. From invariance un- 
der space-time translations, the BS wave function has the 
general form, 

d 4 p 



Xaf}{x,y;P) 



iPX 



(2tt)< 



-W*-v) Xafl (p,P), (A.4) 



with X = (1 — a)x + ay, p := ap a — (1 — a)pp, and 
P = Pa + P/3, where p a , pp denote the momenta of the 
outgoing quarks in the Fou rier transform Xotp{Pot->Pp) P) 
of Xap(x, y; P), see Fig. A.l. One thus has the relation, 



■e-^-^x^P), (A.6) 



X*p(p,P) ■=Xa0(p+(l-cr)P,-p + aP;-P)\ Po=up .(A.5) 

In the definition of the conjugate amplitude, the conven- 
tion 

with X = <jx+ (l—a)y, ensures that hcrmitian conjugation 
from Eq. (|A.2[ ) yields, 

Xa p(P,P) = (7oX t (p,P)7o) a/3 • (A.7) 

In the conjugate amplitude XafiiPi P)t the definition of 
relative and total momenta corresponds to p = (1— <j)p' a — 
op'p and P = —p' a — p'fj for the outgoing quark momenta 



X a i3(p' a ,p'f3;P) = (7oX t ("P^-K; p )7o) Q 







(A. 



c./., Fig. A.l. Note here that hermitian conjugation im- 



plies for the momenta of the two respective quark legs, 
p a — > — p'p, and pp — > — p' a , which is equivalent to a <-> 
(1 — a) and P — > —P. Besides the hermitian conjuga- 



tion of Eq. (A.7), one has from the antisymmetry of the 
wave function, Xup(x,y;~P) = —x/3a(y,x;P). For the cor- 
responding functions of the relative coordinates/momenta, 
this entails that a and (1 — a) have to be interchanged in 
exchanging the quark fields, 

X (x,P)= -X T (-x,P)U (1 _ CT) , (A.9) 

x(p,p) = -x T (~P,P)\ a ^ a) ■ 

This interchange of the momentum partitioning can be 
undone by a charge conjugation, from which the following 
identity is obtained, 

X T ( P , P) = -Cx(-P, -P)^ 1 . (A.10) 
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Fig. A.l. Momentum definitions in diquark amplitudes. 



This last identity is useful for relating x to % in Euclidean 
space. In particular, this avoids the somewhat am bigu ous 
definition of the conjugation following from Eq (A.2) in 
Euclidean space with complex bound-state momenta. 

One last definition for diquark amplitudes concerns the 
truncation of the propagators S attached to the quark legs, 
thus defining the amputated amplitudes x an d X by 



XcfifaP) = (S( Pa )x(p,P)S T (p p )) a0 

x«p(p,p) = {s T (-p' a m P ,p)s(-p' )) 



a/3 



(A.ll) 
(A.12) 



With the definitions above, the same relations hold for the 
amputated amplitudes, in particular, 



x (p,P) = joxHp,Pho 

x(p,p) = -x t (-p,p)\ 



(A.13) 



>(!-*)■ 



In Sec. H the antisymmetric Green function G^ for the 
disconnected propagation of identical quarks with propa- 
gator S(p), 



S a0 (aP-p)S^ 5 ((l-a)P + p) 



(A.14) 
crossed term , 



was used to derive the normalization condition for the 
diquark amplitudes. This notation is somewhat sloppy. In 
particular in the second term proportional to S 4 (p+q), rep- 
resenting the crossed propagation with exchange of the ex- 
ternal quark lines, one may use either p or —q in the argu- 
ments of the propagators. Exchanging one for the other is 
possible only with, at the same time, exchanging a «-> 1— a 
as well, however. Momentum conservation entails that in- 
coming and outgoing quark-pairs in successive correlation 
functions can only be connected if, besides the relative 
momenta p = ±q, also the momentum partitioning vari- 
ables of the pairs match. We can include this condition 
explicitly by introducing temporarily the dimensionless a' 
and a (both in [0, 1]) for the incoming and outgoing pair 
respectively, thus writing, 

G ( ai,0s(P' CT > **, P) = (27r)V(p - q) 5{a - a') (A.15) 
S af} {aP-p)S^((l-(T)P+p) 
-{2Tr) A 5 i {p + q)5{a + a l -I) 

S aS (aP-p)S^{{l-a)P+p) . 



The inverse G^ 1 can then be defined as, 

Gi%s(P, <1, <?', p )=\{ (^Hp ~ (?) 5{a - a') (A.16) 
S-i(aP-p)S^\(l-a)P + p) 
-(2ir) 4 6 4 (p + q)S(a + a'-l) 

S- s 1 (aP~p)S^((l-a)P + p)) , 

giving the exchange antisymmetric unity in the space of 
(identical) 2-quark correlations upon (left or right) multi- 
plication with G' ' , 
l 

/ J d ° Gi °ipJP>^ k >°>nG$]phk^,q,<?',P) = 

o 

i (Sap S lS (2vr) 4 J 4 (p - q) S(a - a') (A.17) 



■SasSrf {2n) 4 d 4 {p + q)6(a 



1 



This is the way the multiplication of 2-quark correlation 
functions for identical quark s used in Sec. B is understood 
properly. Either Eq. ( A.15 ) or Eq. ( A. 16| ) can be used 
in the derivation of the normalization conditi on for the 
diquark amplitude x f° r identical quarks, Eq. (|2.9|). 



B Supplements on the Nucleon BSE 

In this appendix we would first like to explore the pos- 
sibility of having exchange-symmetric arguments in the 
diquark amplitudes of the quark-exchange kernel of the 
nucleon BSE. Consider the invariant x\ for the relative 



momentum in the incoming diquark in the kernel (3.6) 



corresponding to the defin ition of the s ymm etric argum ent 
of P(x) given in Eq. fl2.12| ). From Eqs. (|J) and ( |3TlD one 
obtains (with er + a = 1 and r\ + f) = 1 determining the 
momentum partitioning within the diquark and nucleon 
respectively), 

Xi = —p 2 — oak 2 — pk + (1 — 3r])Pp+ {2aafj—r])Pk 

+(r7(l-2?7) - aaf)P 2 . (B.l) 

On e veri fies readily that for a = (1 — 2?7)/ (1 — 77), according 
to (3.13), the prefactor of the term oc P 2 vanishes, and 
that of the term oc Pk becomes a(l— 3rj) leaving only the 
choice 77 = 1/3 for a P-independent x±. One might now 
argue that, from the antisymmetry considerations alone, 
we should be free to add an arbitrary term proportional 
to the square of the total diquark momentum in the def- 
inition of x. We will show now that such a redefinition 
cannot lead to a P-independent x either (for values of 
77 different from 1/3). Here, the diquark momentum is 
Pd = q + Pa = f]P ~ k and with 



x\ := x\ 

2 D 2 



CP 2 



Pi = fj z P z - 2rjPk 
one finds that in order to have no terms oc P 2 in x\, 



(B.2) 



C= -ff 2 (77(1 — 277) -ctctt) 2 ) 



(B.3) 
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With this C, however, one has 

Xi = -p 2 - r]7j- 2 {l-2ri)k 2 - pk + (l-3r]){Pp + tit)- 1 Pk) , 

independent of a. This shows that the symmetric argu- 
ments of the diquark amplitudes x and, analogously, \ 
can quite generally be independent of the total nucleon 
momentum P only for r\ — 1/3. 

The remainder of this appendix describes the struc- 
ture of the nucleon BSE for the bound-state of scalar di- 
quark and quark in some more detail. The form of the 
bou nd-s tate pole in the scalar-fermion 4-point function, 
Eq. (3.1), implies that the corresponding bound-state am- 
plitudes obey, 



^(p,P n )A+(P n 
A+{P n )$(p,P n 



1/>(P, Pn) , 



(B.4) 



with A + (P n ) = (rf n + M n )/2M n . Therefore, the ampli- 
tudes can be decomposed as follows: 



Si(p,P n )A + (P n ) + S 2 (p,P n )S(p,P n )A + (P n ) , 

$(P,Pn) = (B.6) 

Si(-p,P n )yl + (P n ) + S 2 (-p,P n )A+(P n )S(-p,P n ) , 

with S(p, P n ) = (^>—pP n /M n )/M n . This simply separates 
positive from negative energy components of the ampli- 
tudes, 



If n A+{P n ) = M n A+(P n ) , 
P n S(p,P n )A+(P n ) = ~M n S(p,P n )A+(P n ) , 
and thus, A+{P n ) 5(p, P n ) A+{P n ) = . 

One furthermore has, 

A+(P n ) S(p, P n )S(p, P n ) A+(P n ) = 



(B.7) 



P 

Ml 



Ml 



(B.8 

A+(P) , 



which allows to rewrite the homogeneous BSE (3.4) in 
terms of 2-vectors S T (p, P n ) := (Si(p,P n ), S 2 (p, Pn)), us- 
ing the kernel (3.6), 



S( P ,Pn) 
1 



2N*J (2tt)« 
with k s = (l—r))P n 

T(p,k,P n ) = \[\ 



(B.9) 

P(x 1 )P(x 2 ) D(k s ) T(p, k, P n ) S(k, P n ) , 



-k 





r\ P n and 



p 

M 2 



Mi 



(B.10) 



tv{S(q)S(k q )A+(P n )} 
tv{S(q)S(k q )A+(P n )E( P ,P n )} 

tl {S(q)S(k q )S(k,P n )A+(P n )} 
tr {S(q)S(k q )S(k, P n )A+(P n )E(p, P n )} 

After performing these traces the transfer to Euclidean 
metric introducing 4-dimensional polar variables is done 



according to the prescriptions, 
p 2 , k 2 -> -p 2 , -k 2 , Pi -> M 2 
pP n -> iM n p y , kP, 
pk — > —kpu(x,y,z) , 

u(x,y,z) 



n ' 

iM„ k z . 



(B-ll) 



yz + x yl — y 2 \/\ — z 2 . 

In these variables, the nucleon BS-amplitudes are func- 
tions of the modulus of the relative momentum and its 
azimuthal angle with the total momentum. The matrix T 
in the kernel, in addition to the moduli p, k and azimuthal 
angles y, z of both relative momenta, also depends on the 
angle u between them, 

S(p,P n ) -> S(p,y) , 

T(p, k, P n ) -> T(p, y, k, z, u(x, y, z)) . (B.12) 

The azimuthal dependence of the amplitudes is taken into 
account by means of a Chebyshev expansion to order TV, 
see, e.g., Ref. J60[, 



N-l 



(B.5) S (p,y) ~ ^Hr5„(p)T„(y), 

n=Q 

2 N 

Snip) = i n ^Sip, Vk)T n {yk) , 

k=l 



(B.13) 
(B.14) 



where the y k — cos 



7r(fe-l/2) 
TV 



are the zeros of the Chebyshev polynomial of degree N. 
Here, Chebyshev polynomials of the 1st kind are used 
with, for later convenience, a somewhat non-standard nor- 
malization To := 1/V2. An explicit factor (— i) n was intro- 
duced in order to obtain real Chebyshev moments S n (p) 
for all n. Analogous formulae are obtained for expansions 
in Chebyshev polynomials of the 2nd kind, which are used 
in Refs. |36l|37||. The nucleon BSE now reads, 



s m {p) = 



with 



(B.15) 



k 3 dk 



N-l 



2N1 / (4tt) 2 ^ 

v ; n=0 



i m n T mn (p, k) S n {k) 



(B.16) 



T mn (p,k) = — J \/l - z 2 dz 
-1 

1 r 2 N ( 

T(p,y k ,k,z,u(x,y k ,z)) T m (y k )T n {z) , 

where x s = k 2 + 2i(l—T])M n kz—(l—ri) 2 M 2 is the invariant 
momentum of the free scalar propagator D of mass m s . 

C Supplements on Nucleon Charge 
Conservation 

The missing step in the explicit verification of the correct 
charges Qp = 1 for the proton and Qn = for the neutron 
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J tp M SMr£ M SMDMlpM 



Wick rotation 



Jj E s E r» E s E D E 4, E 

+ residue terms 



Fig. D.l. Interrelation of matrix elements in Minkowski and Euclidean space. The i nteg ral sign is 
shorthand for the four-dimensional integration over the relative momentum k, see Eq. (D.l). 



(C.l) 



is to prove that 

N q — N D = 2N X + 3iV P 

which is equivalent to Qn — 0, see Eqs. ( |6.27|) to ( |6.29 ) 
in Sec. [| To this end, note that from Eqs. ( ^4.5| ) and ([4.6) 
the l.h.s above can be written as, 

^ (C2) 



N, 



N °=2Mj 
pv 



tr 



tr 



2M n 



(27T) 4 

d 
d 4 k 



S- l (k q )D~ l {k a ) Tjj(k,P) 



(2tt)* 
d 



S-^kJD-^kJMktP) 



^(-k,P)S- L (k q )D- L (k s ) i>(k,P) 



total derivative 



with k q = r)P + k and k s = fjP — k as in the previous 
sections. A surface term which vanishes for normalizable 
BS wave functions i\) and ip was not given explicitly. Using 

the BSEs for if = S^D- 1 ^ and % = ^S^D' 1 , c.f., 
Eqs. (5.7) and (J5^), in analogy to the proof of current 
conservation in Sec. |)[ one obtains, 



tr 



d 4 p d 4 k 



2M n J (2tt) 4 (2tt) 



(C.3) 



From Eq. ( p.6[ ) for the explicit form of the exchange kernel 
it follows that 

d : _d_\ „, , _ 1 



dk^ J 



K(p,k,P) = - 



2m 



(C.4) 



2 ( q^ s ^ ) p (-pI)p(-pI) 



S(q) 



^P[{-p\)P{-pl) - 3p 2l ,Pi(-p 2 1 )P'(-p 2 2 ) 



since q= (l-2rj)P-p-k,pi = -{l-3r))P/2+p + k/2 and 
P2 = (l — 3r))P/2 — p/ 2 — k. Co mpa ring to the definitions 
of N x and N P , Eqs. (17) and p^ ) in Sec . [| resp ectiv ely, 
we see that Eq. (C.4) inserted in Eq. ( |C.3|) gives (C.l) as 
required. 



D Calculation of the Impulse Approximation 
Diagrams 

In this appendix we discuss the difficulties in the formal 
transition from the Minkowski to the Euclidean metric 
which are encountered in the connection between Bethe- 
Salpeter amplitudes tp and Bethe-Salpeter wave functions 
ip in a general (boosted) frame of reference of the nucleon 
bound state. 

As the generic exampl e fo r this discussion, we choose 
the first diagram in Figure |6~T| which describes the impulse- 
approximate contribution arising from the coupling of the 
photon to the quark within the nucleon, (J^ 1 ) according 
to Eq. (|6.10| ) with Eq. ( |6.11 ). Please refer to Fig. 6.1 and 
Table |6.15| for the momentum definitions employed herein. 

In the Mandelstam formalism, such matrix elements 
between bound states are related to the corresponding BS 
wave functions in the in Minkowski space, here to the nu- 
cleon BS wave functions V'a/ -0 Upon the transition to the 
Euclidean metric, the corresponding contribution to the 
observable, here to the nucleon form factors, is determined 
by the "Euclidean BS wave function" tp^.hi the rest frame 
of the nucleon bound state this transfer from M — ► E of 
the BS wave functions commutes with the replacement of 
the wave functions by the truncated BS amplitudes; that 
is, unique results are obtained from the Euclidean contri- 
butions based on either employing the Minkowski space 
wave functions or the BS amplitudes which are related by 
the truncation of the propagators of the constituent legs, 
here V>m = S^D^ipM, or, vice versa, ipM = SmD m ^m- 
At finite momentum transfer Q 2 one needs to employ 
BS wave functions in a more general frame of reference, 
here we use the Breit frame in which neither the incoming 
nor the outgoing nucleon are at rest. As described in Sec. 

the "Euclidean wave function" ipE in this frame is ob- 
tained from the solution to the BSE in the rest frame by 
analytic continuation, in particular, by inserting complex 
values for the argument of the Chebyshev polynomials, 
see Eqs. ( |6.33 ). This corresponds to the transition from 
left to righ t indicated by the arrow of the upper line in 



Figure D.l. 

In the analogous transition on the other hand, when 
the truncated BS amplitudes are employed, the possible 
presence of singularities in the propagators of legs has to 

5 In the following the subscript m stands for definitions in 
Minkowski space and E for the corresponding ones in Euclidean 
space. 
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be taken into account explicitly. In the present example, 
these are the simple particle poles of the propagators of 
the constituent quark and diquark that might be encircled 
by the closed path in the fco-integration. The correspond- 
ing residues have to be included in the transition to the 
Euclidean metric in this case, which is indicated in the 



lower line of Fig. D.l 



The conclusion is therefore that the naive relation be- 
tween BS amplitudes and wave functions can not be main- 
tained in the Chebyshev expansion of the Euclidean spher- 
ical momentum coordinates when singularities are encoun- 
tered in the truncation of the legs. Resorting to the Min- 
kowski space definitions of BS amplitudes vs. wave func- 
tions, however, unique results are obtained from either em- 
ploying the domain of holomorphy of the BS wave func- 
tions in the continuation to the Euclidean metric (with 
complex momenta) or, alternatively and technically more 
involved, from keeping track of the singularities that can 
occur in the Wick rotation when the truncated amplitudes 
and explicit constituent propagators are employed. 

The rest of this section concerns the description of how 
to account for these singularities which, for our present 
calculations employing constituent poles for quark and di- 
quark, give rise to r esid ue terms as indicated in the lower 
right corner of Fig. 



D.l 



To this end consider the quark contribution to the ma- 
trix eleme nts of the elect romagnetic current which, from 
Eq. (|6T0|) with Eq. ( PH ), is given by 



d 4 k 



(D.l) 



(2tt)4 

- f)Q, p') D(k s )S{ Pq )r» uark s{k q )4>{k, p) . 

We are interested in the location of the propagator poles 
herein. For these poles, solving the corresponding quadratic 
equation for the zer oth c omponent of th e relative momen- 
tum k , from Table 6.15 and Eqs. (6.3C), yields 



kg - m\ - ie = 
k° 

"'pole,! 
Pa ~ m a ~ * e = <^ 



-77 loq ± W{m q ,(k- v /2Cl) 2 



(D.2) 



(D.3) 



k° P oie,2 = -WQ ± VF(m g ,(k-(, ? /2-l)Q) 2 ) 

(D.4) 



ki - mi - ie = <^ 

1,0 

% pole,3 



fcL., = f,uj Q ± ^(m s ,(k + r)/2Q) 2 ) 



with W(m, k 2 ) = Vk 2 + m? — ie (and r] + f]= 1). 

For Q = 0, i.e., in the rest frame of the nucleon in 
which the Bethe-Salpeter equation was solved, the naive 
Wick rotation is justified for 1 — ^ < rj < ^rj-, since there 
always is a finite gap between the cuts contained in the 
hypersurface Refc°=0 of the Refc - Imfc - k space. As Q 
increases, these cuts begin are shifte d alo ng both, the k 3 
and the fc°-axis, as sketched in Figure D.2. This eventually 



amounts to the effect that one of the two cuts arising from 
each propagator crosses the Im fc°-axis. As indicated in 
the figure, Wick rotation k° — > ik A is no longer possible 
for arbitrary values of k 3 without encircling singularities. 



Cut *, 



encircled singularities 



. Re k° 




Fig. D.2. Location of the relevant singularities in the impulse- 
approximate quark contribution to the form factors. The rel- 
ative momentum k is the in tegra tion variable in the loop dia- 
gram corresponding to Eq. (D.l). 



The corresponding residues thus lead to 
d A k 



(J?) 



1q 



(2tt)4 



(D.5) 



+i 



i>{-k - f,Q, P 1 ) D(k s )S( Pq )r»S(k q )^(k, P) 

d 3 k 



k A 



(2tt) 3 

r£Jies(S(k q ))i;(k* oleil ,k,P) 
+ analogous terms for S(p q 



k-f,Q,P')D(k s )S( Pq ) 



and D(k s 



upon transforming Eq. (D.l) to the Euclidean metric. Here, 
the residue integral is evaluated at the position of the pole 
in the incoming quark propagator S(k q ) on the Euclidean 
ft 4 -axis 



Koie,i = -*vVM 2 + Q 2 /4 + tW(m q , (k - V /2 Q) 2 ) , 

where Res(S(k q )) denotes the corresponding residue, and 
the abbreviation 



ynuQ 



^(m 9 ,(k-77/2Q) 2 )) 



was adopted to determine the integr ation domain for which 
the encircled singularities of Figure [D.2| contribute. 

Analogous integrals over the spatial components of the 
relative momentum k arise from the residues correspond- 
ing to the poles in the outgoing quark propagator S(p a ) 
and the d iquark propagator D(k s ) as given in Eqs. (D.2) 
and |51|). 

One verifies that the se cuts (as represented by the 
shaded areas in Fig. p.2j ) never overlap. Pinching of the 
deformed contour does not occur, since there are no anoma- 
lous thresholds for spacelike momentum transfer Q 2 in 
these diagrams. 
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Proton Electric Form Factor 

m =0.58 GeV, m d =0.66 GeV, P(p): dipole 
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Fig. D.3. The impulse-app roxia mte contribution, correspond- 
ing to the diagrams of Fig. 



6.1 



to the electric form factor of 
the proton (employing the dipole diquark amplitude). Results 
of the BS wave function calculations for the 77- values 1/3 and 
0.4 are compared to the respective amplitude plus residue cal- 
culations. The results are normalized to the latter with 77=1/3. 



In Figure D.3 we compare results for the electric form 
factor of the proton employing the nucleon BS amplitudes, 
together with the procedure to account for the singulari- 
ties as outlined above, with the corresponding wave func- 
tion calculations. The n — 2 dipole diquark amplitude of 
Sec. ^ is employed herein once more. For r\ = 1/3 the 
Chebyshev expansion of the BS wave function to 9 orders 
still turns out insufficient to provide for stable numerical 
results. This is due to being too close to the range in 77 
that requires proper treatment of the diquark pole contri- 
bution to the nucleon BSE which, with the present value of 



the mass m s /M = 0.7 is the case for r\ < 1 — m s /M = 0.3. 
The considerably weaker suppression of higher orders in 
the Chebyshev expansion of the BS wave function as com- 
pared to the expansion of the BS amplitude enhances the 
residual ?y-dependence of the observables obtained from 
the former expansion at a given order, in particular, when 
it has to reproduce close-by pole contributions in the con- 
stituent propagators. The impulse-approximate contribu- 
tions to Ge deviate substantially from those employing 
the BS amplitude and residue calculations in this case. On 
the other hand, for values of the momentum partitioning 
variable which are a little larger than 1/3 such as r\ = 0.4 
used in the other results of Fig. D.3, unique results are 
obtained from both procedures. Both, the BS wave func- 
tion and amplitude calculation are in perfect agreement 
for values of the momentum partitioning that are closer 
to the middle of the range allowed to r\. 

For the seagull and exchange q uark contributions cor- 
responding to the diagrams of Fig. 3.2 an analogous anal- 
ysis of the singularity structure is considerably more com- 
plicated. The explicit inclusion of their contributions which 
allowed the calculation based on the BS amplitude expan- 
sion also of these diagrams is numerically too involved. For 
these contributions to the form factors we have to resort 
to the BS wave function calculations. Unlike the impulse- 
approximate contributions we find, however, that the de- 
viations in the results for the exchange quark and seagull 
contributions for 77 = 1/3 and r\ — 0.4 are smaller than the 
numerical accuracy of the calculations and thus negligible. 
We attribute this to the fact that the exchange quark and 
seagull contributions to the form factors tend to fall off 
considerably faster with increasing Q 2 than those of t he 



impulse approximation. This can be seen in Figure 6.6 



Small residual 77-dependences are observed for the mo- 
mentum transfers above 3GeV 2 also in the otherwise sta- 
ble calculations. These give rise to deviations in the results 
for the form factors of at most 4% which decrease rapidly 
and become negligible at lower Q 2 . 
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